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We calculate the dynamic structure factor S(u,q) of spinless fermions in one dimension with 
quadratic energy dispersion k 2 /2m and long range density-density interaction whose Fourier trans- 
form f q is dominated by small momentum-transfers q < qo <C &f- Here qo is a momentum-transfer 
cutoff and kF is the Fermi momentum. Using functional bosonization and the known properties of 
symmetrized closed fermion loops, we obtain an expansion of the inverse irreducible polarization to 
second order in the small parameter qo/^F- In contrast to perturbation theory based on conven- 
tional bosonization, our functional bosonization approach is not plagued by mass-shell singularities. 
For interactions which can be expanded as /„ = f + /oY/2 + 0(q 4 ) with /q ^ we show that 
the momentum scale q c — l/]m/o'| separates two regimes characterized by a different (/-dependence 
of the width 7 q of the collective zero sound mode and other features of S(lo, q). For q c <C q <C kF 
all integrations in our functional bosonization result for S(ui, q) can be evaluated analytically; we 
find that the line-shape in this regime is non-Lorentzian with an overall width 7 9 oc q 3 /(mq c ) and 
a threshold singularity [(uj — ln 2 (cj — cj q ')]^ 1 at the lower edge uj — > uj^ — vq — jg, where v is 
the velocity of the zero sound mode. Assuming that higher orders in perturbation theory transform 
the logarithmic singularity into an algebraic one, we find for the corresponding threshold exponent 
fiq — 1 — 2rj q with r) q oc q 2 jq 2 . Although for q < q c we have not succeeded to explicitly evaluate our 
functional bosonization result for S(u,q), we argue that for any one-dimensional model belonging 
to the Luttinger liquid universality class the width of the zero sound mode scales as q 2 /m for q — > 0. 

PACS numbers: 71.10.Pm, 71.10.-w 



I. INTRODUCTION 

Recently several authors have calculated the dy- 
namic structure factor S(u>, q) in the Luttinger liq- 
uid phase of model systems for interacting fermions 
with non-linear energy dispersion in one spatial dimcn- 
sion i 1 2 i 3 1 4 i 5 i 6 i 7 i 8 i 9 i io^i i 2 i 3 i 4 i 5_ Mathematically, S{u>,q) 

is defined as the spectral density of the density-density 
correlation function, 

S(u, q) = Jdtjdx e*^-*") (5p(x, i)Sp(0, 0)) , 

(1.1) 

where Sp(x, t) is the operator representing the deviation 
of the density from its average. The dynamic structure 
factor can be directly measured via scattering experi- 
ments probing density-density correlations of the system. 
It is therefore important to have quantitatively accurate 
theoretical predictions for the line-shape of S(tv, q). 

Although there is general agreement that in the 
Luttinger liquid regime of one-dimensional interacting 
fermions S(u,q) exhibits for small frequencies lj and 
wave-vectors q a narrow peak associated with the col- 
lective zero sound (ZS) mod o 16 : 17 , a quantitative under- 
standing of the precise line-shape of the ZS resonance in 
generic non-integrable models is still lacking. The spec- 
tral line-shape is expected to depend on non-universal 
parameters of the model under consideration, such as the 
non-linear terms in the expansion of the energy dispersion 



efc around the Fermi momentum kp, or the coefficients 
in the expansion of the Fourier transform f q of the inter- 
action for small momentum-transfers q. Because these 
parameters correspond to couplings which are irrelevant 
(in the renormalization group sense) at the Luttinger liq- 
uid fixed point, the line-shape of S(oj,q) cannot be ob- 
tained using standard field-theoretical methods, such as 
field-theoretical bosonization, which has otherwise been 
very successful to obtain the infrared properties of Lut- 
tinger liquid o 18 ' 19 : 20 ' 21 . Recall that the crucial step in 
the bosonization approach is the linearization of the en- 
ergy dispersion around the Fermi points, €k F +q — £fe F ~ 
VFq, where vf is the Fermi velocity. If in addition the 
Fourier transform f q of the interaction is non-zero only 
for momentum-transfers q <C kp, we arrive at the ex- 
actly solvable Tomonaga-Luttinger model (TLM), whose 
bosonized hamiltonian is non-intcractin g 18 i 19 i 20 i 21 . As a 
consequence, the dynamic structure factor of the TLM 
has only a single 5-function peak corresponding to a 
collective ZS mode with infinite lifetime. For spinless 
fermions with long-range density-density interaction f q 
one obtains for small q, 

Stj,u{^,q) = Z q S(uj - v \q\) , (1.2) 

where the velocity vo and the weight Z q of the collective 
ZS mode can be written as 



vo/vf = y/l + go, (1-3) 
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= v F q 2 = M 
9 27rt;o |g| 27i7T+Sb' V ' 7 

For later convenience we have introduced the relevant di- 
mensionless interaction at vanishing momentum-transfer, 

go = t'o/o, (1-5) 

where uq = l/(nvp) is the non- interacting density of 
states at the Fermi energy. 

The question is now how the line-shape of S(uj,q) 
changes if wc do not linearize the energy dispersion. 
There have been many recent attempts to find an an- 
swer to this question. Roughly, the proposed methods 
can be divided into four different categories: 

1. Conventional bosonization. The established ma- 
chinery of conventional bosonizatio n 18 i 19 i 20 ' 21 has been 
used in Refs. [llMl-ij to calculate the dynamic struc- 
ture factor of Luttinger liquids. Expanding the en- 
ergy dispersion around k = kp beyond linear order, 
tkp+q ~ ^kp+vpq+q 2 / {2m), the quadratic term q 2 /(2m) 
gives rise to cubic interaction vertices proportional to 
1 /m in the bosonized model 2 ^. Hence, bosonization maps 
the original (unsolvable) fermionic many-body problem 
onto another unsolvable problem involving bosonic de- 
grees of freedom. The hope is that perturbation theory 
for the effective boson model is well-defined and more 
convenient to carry out in practice than in the original 
fermion models. Unfortunately, this strategy fails for 
the calculation of S(ui, q), because already to second or- 
der in 1/m one encounters singular terms proportional 
to l/(u;±voq), which become arbitrarily large as the fre- 
quency approaches the mass-shell uj — ► ±voq. Some time 
ago Samokhini proposed a simple regularization proce- 
dure of these mass-shell singularities which we shall re- 
view in Sec. El Assuming a Lorentzian line-shape, he 
found that for q — > most of the spectral weight is 
smeared out over an interval of width q 2 /m. Although 
this estimate for the width of the ZS resonance was later 
confirmed by various other calculations ^ 9 ' 11 ^ 11 ' 14 , the as- 
sumption of a Lorentzian line-shape turns out to be in- 
correct. It would certainly be more desirable to have a 
controlled method of re-summing the interaction in the 
bosonized hamiltonian to infinite orders such that the un- 
physical mass-shell singularities arc properly regularized; 
apparently this problem has not been solved so far. We 
shall further elaborate on these mass-shell singularities 
in Sees. Eland US 

2. Re-summing fermionic perturbation theory via an 
effective hamiltonian. Because of the above mentioned 
problems inherent in standard bosonization, it seems bet- 
ter to set up the perturbation expansion in terms of the 
original fermionic degrees of freedom using diagrammatic 
techniques. In this approach, it is convenient to first cal- 
culate the polarization function T\(iu, q) for imaginary 
frequencies and then use the fluctuation-dissipation the- 
orem to obtain the dynamic structure factor, 

S{uj, q) = 7r _1 Imn(o; + iO, q) . (1.6) 



For simplicity, we shall focus on the limit of vanishing 
temperature throughout this work. For long-range inter- 
actions whose Fourier transforms f q are dominated by 
small wave- vectors q <C kp, one usually avoids the di- 
rect expansion IL(u>, q) in powers of the bare interaction, 
but instead expands its irreducible part II* (u>, q) which 
is defined via 

n- 1 (cu,q) = f q + n- 1 (uj,q). (1.7) 

In a recent paper, Pustilnik et alA did not follow this 
standard approach, but expanded the full (i.e., reducible) 
polarization IL(u,q) in powers of the bare interaction. 
They found already at the first order in the bare in- 
teraction that the correction to S(uj,q) diverges loga- 
rithmically if w approaches a certain threshold edge u)~ 
from above. The authors of Ref. [i| then proposed a re- 
summation procedure of the most singular terms in the 
perturbation series to all orders using an effective hamil- 
tonian constructed in analogy with the X-ray problem. 
In this way, they succeeded to transform the logarith- 
mic threshold singularity into an algebraic one, char- 
acterized by a certain momentum-dependent threshold 
exponent. The spectral line-shape can therefore not be 
approximated by a Lorentzian as implicitly assumed by 
Samokhin 1 ; on the other hand, Samokhin's result that 
the overall width of the ZS resonance scales as q 2 /m was 
confirmed by Ref. 0]. However, Pustilnik et alX did not 
explicitly analyze the higher-order terms in the perturba- 
tion series to demonstrate that the logarithmic singular- 
ity encountered at the first order can really be re-summcd 
to all orders to yield an algebraic singularity. Moreover, 
they did not keep track of the (finite) renormalization 
of the ZS velocity v, which determines the precise en- 
ergy scale of the collective ZS resonance and its position 
relative to the energy of the single-pair particle-hole con- 
tinuum, which a priori need not be identical. 

3. Integrable models. Another method to calculate the 
dynamic structure factor of Luttinger liquids is based 
on the analysis of exactly solvable models belonging 
to the Luttinger liquid universality class, such as the 
XXZ-chain 9 - ' 10 ' 11 or the Calogero-Suthcrland models. 
These calculations have confirmed the results obtained 
by Pustilnik et alA for generic (not necessarily inte- 
grable) one-dimensional Luttinger liquids: The spectral 
line-shape is non-Lorentzian, exhibits algebraic threshold 
singularities, and the weight is smeared over a frequency 
interval proportional to q 2 jm for q — > 0. But since the 
re-summation procedure of Ref. [J] is not rigorous and 
higher-order terms in the perturbation series have not 
been explicitly analyzed, one cannot exclude the possi- 
bility that the algebraic threshold singularities are a spe- 
cial feature of integrable models, and that in generic non- 
integrable models the higher order terms in the perturba- 
tion series do not conspire to transform logarithmic sin- 
gularities into algebraic ones. Note also that the effective 
two-body interaction in the spinless fermion model ob- 
tained from the XXZ-chain via the usual Jordan- Wigncr 
transformation involves also momentum-transfers of the 



order of kp. This model is therefore different from the 
forward scattering model with quadratic dispersion con- 
sidered here, where the Fourier transform of the density- 
density interaction f q is only finite for q <C kp. Appar- 
ently, an exactly solvable model with non-linear energy 
dispersion and density-density interaction f q involving 
only small momentum transfers and / g =o > does not 
exist. Although the momentum dependence of the inter- 
action is irrelevant in the renormalization group sense, 
the line-shape of S(u),q) is essentially determined by ir- 
relevant couplings, so that models with different sets of 
irrelevant couplings might also exhibit different spectral 
line-shapes. 

4- Functional bosonization. This is an alternative 
method of describing fcrmionic many-body systems with 
dominant forward scattering in terms of bosonic degrees 
of freedom. In the context of the TLM, the functional 
bosonization idea has been introduced by Fogedby 2 ^ and 
by Lee and Cher* 2 ^. Later this technique has been used 
to bosonize interacting fermions with dominant forward 
scattering in arbitrary dimensions 2 ^, and to estimate the 
effect of the non-linear energy dispersion on the single- 
particle Green function^. For a review of this approach 
see Ref. [27j, where the advantages of this method for 
calculating the dynamic structure factor have already 
been advocated. Like in conventional bosonization, in the 
functional bosonization approach the non-linear terms in 
the energy dispersion give rise to interaction vertices in 
the effective bosonized action of the system. However, 
the interaction vertices in functional bosonization arc 
rather different from the vertices due to the non-linear 
dispersion in conventional bosonization. In fact, the in- 
teraction vertices in functional bosonization can be iden- 
tified diagrammatically with symmetrized closed fermion 
loops, which can be calculated exactly for quadratic dis- 
persion in one dimensio n 28 i 29 i 30 . While in conventional 
bosonization a quadratic energy dispersion gives rise to 
cubic vertices in the bosonized Hamiltonia n 19 ' 22 , within 
functional bosonization a quadratic dispersion leads to 
infinitely many vertices involving an arbitrary number 
of boson fields. The fact that perturbation theory for 
S(u>, q) based on functional bosonization is different from 
perturbation theory based on conventional bosonization 
is obvious if one considers the non- interacting limit: while 
functional bosonization yields the exact free polarization 
Ho(u),q), conventional bosonization produces an expan- 
sion of IIo(w, q) in powers of 1/m, which in practice has 
to be truncated at some low order, leading to unphysical 
mass-shell singularities. 

In Ref. [IH two of us have used the functional bosoniza- 
tion approach to calculate the width j q of the ZS mode 
in a generalized Tomonaga model with quadratic energy 
dispersion. To estimate the effect of non-linear energy 
dispersion on the dynamic structure factor, we have trun- 
cated the expansion of the inverse irreducible polariza- 
tion at the first order in an expansion in powers of the 
Gaussian propagator of the boson fields, which can be 
identified with the effective screened interaction within 
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FIG. 1: Diagrammatic definition of the screened interaction 
within random phase approximation. The thin wavy line de- 
notes the bare interaction and the solid arrows represent non- 
interacting fermionic single-particle Green functions. 




(a) (b) (c) 

FIG. 2: Corrections to the irreducible polarization in an ex- 
pansion to first order in powers of the RPA interaction. 



random phase approximation (RPA) defined in Fig. [T] 
To this order, the simple first-order Hartree contribu- 
tion to the bosonic self-energy in the functional bosoniza- 
tion approach (the corresponding Feynman diagram is 
shown in Fig. \E\ (a) in Sec. IIV|) is in fermionic language 
equivalent to the sum of the three first-order interac- 
tion corrections to the irreducible polarization shown in 
Fig. Functional bosonization thus consistently sums 
self-energy corrections (diagrams (a) and (b) in Fig. ^ 
and vertex corrections (diagram (c) in Fig. [2]) of the 
underlying fermion problem. Actually, the interpreta- 
tion of the inverse irreducible polarization as the self- 
energy of the effective boson theory obtained via func- 
tional bosonization suggests that one should always ex- 
pand the inverse II^ 1 (ui, q) in powers of the relevant small 
paramete r 12 i 27 . 

Unfortunately, it is not consistent to truncate the ex- 
pansion of n~ 1 (w, q) at the first order in the RPA inter- 
action, so that the result j q oc q 3 for the ZS damping 
obtained in Ref. [l2| cannot be trusted. We shall ex- 
plain this in more detail in Sec. IIV1 where we construct 
a systematic expansion of IT^ 1 ^, q) in powers of bosonic 
loops using functional bosonization and show that for 
our forward scattering model there is a large interme- 
diate regime q c < q -C kp where indeed 7 g oc q 3 /(mq c ). 
The momentum scale q c is determined by the momentum 
dependence of the interaction f q , see Eq. (| 1 . 16[) below. 
Due to the complexity of the integrations, in the regime 
q <C q c we have not been able to evaluate our functional 
bosonization result for S(u,q). However, at q w q c our 
expression for j q matches the result 7 g oc q 2 /m obtained 
by several other authors for different model systems for 
Luttinger liquid*^^. We therefore believe that quite 
generally for any model belonging to the Luttinger liquid 
universality class the width of the ZS resonance asymp- 
totically scales as q 2 for q — > 0. 
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Let us now briefly define our model and introduce 
some useful notation. We consider non-relativistic spin- 
less fermions interacting with long-range density-density 
forces in one spatial dimension. The Euclidean action is 



S[C, C] = S [C, C} + - I fqP-QPQ ■ 

1 Jo 



(1.8) 



where the non-interacting part can be written in terms 
of Grassmann fields ck and ck representing the spinlcss 
fermions as follows, 

S [c,c] = - / (iu> - e k + [i)c K c K - (1-9) 
■Ik 

Here, /j, is the chemical potential and the energy disper- 
sion is assumed to be quadratic, 



a- 



2m 



The composite field 



PQ 



K 



crcr+q 



(1.10) 



(1.11) 



represents the Fourier components of the density. The 
collective label K = (iw, k) denotes fermionic Matsub- 
ara frequencies ilu and wave- vectors fc, while Q = (iQ, q) 
depends on bosonic Matsubara frequencies iQ. The cor- 
responding integration symbols are J„ = (f3V)~ 1 ^ w fc , 
and Jq = (PV)^ 1 q' wnere P is the inverse tem- 
perature and V is the volume of the system. Eventu- 
ally, we shall take the limit of infinite volume V — > oo 



dudk 



and 



and zero temperature (3 — > oo, where J K = J 

Jq = J (2-11$ • ^ e assume that the Fourier transform f q 
of the interaction is suppressed for momentum-transfers 
q exceeding a certain cutoff go ^ &f- For explicit calcu- 
lations it is sometimes convenient to use a sharp cutoff 12 , 



/, = fo&(qo - \q\) • 



(1.12) 



However, as will be discussed in detail in Sec.[V] the van- 
ishing of all derivatives of f q at q = eliminates an im- 
portant damping mechanism, so that it is better to work 
with a more realistic smooth cutoff, such as a Lorcntzian, 



fo 



i + q 2 /Qo 



(1.13) 



Throughout this work we assume that the momentum- 
transfer cutoff qo (which for Lorentzian interaction can be 
identified with the Thomas- Fermi screening wave- vector) 
satisfies 



(1.14) 



The precise form of f q is not important for our purpose, 
as long as for small q we may expand 

U = fo + l:foq 2 + 0(q 4 ), with/oVO. (1.15) 



By dimensional analysis, we may use the second deriva- 
tive /o' of the Fourier transform of the interaction to con- 
struct a new momentum scale 



1 



H/o'l ' 



(1.16) 



which will play an important role in this work. Note 
that for Lorentzian cutoff /q = — 2/o/q 2 < and 
q c = qio/(2m/o), but in general the momentum scale q c is 
independent of the ultraviolet cutoff q . We assume that 



q c < qo <C k F . 



(1.17) 



For simplicity, we shall refer to the model defined above 
as the forward scattering model (FSM). If we further sim- 
plify the FSM by linearizing the energy dispersion around 
the two Fermi points, e±k F +q — £k F ~ ±vpq and by ex- 
tending the linear dispersion at each Fermi point to the 
infinite line — oo < q < oo, then the FSM reduces to the 
spinlcss TLM with dimcnsionless forward scattering in- 
teractions g2 = 34 = go in "g-ology" -notation^. In con- 
trast to the TLM, the FSM does not require ultraviolet 
regularization, because the quadratic energy dispersion 
in one dimension renders all loop integrations ultravi- 
olet convergent. Hence the usual problems associated 
with the removal of ultraviolet cutoffs and the associated 
anomalies^ simply do not arise in the FSM. 

To conclude this section, let us give a brief outline 
of the rest of this work. In Sec. |H] we shall discuss 
the dynamic structure factor of the FSM within the 
RPA; although in this approximation the ZS mode is not 
damped, it is still instructive to start from the RPA be- 
cause it allows us to understand the origin of the mass- 
shell singularities encountered in conventional bosoniza- 
tion. In Sec. HIH we outline the functional bosonization 
approach to the FSM, which we then use in Sec. IIVI to 
derive a self-consistency equation for n~ 1 (w, q) which 
does not exhibit any mass-shell singularities. In Sec. [V] 
we present an evaluation of this expression for sharp 
momentum-transfer cutoff (|1.12[) . while in Sec. I VII we 
consider a general interaction f q of the type (jl.lSp . We 
also present explicit results for the spectral line-shape of 
S(u>, q) and the ZS damping. In Scc. lVIll we briefly sum- 
marize our main results and point out some open prob- 
lems. There are four appendices with technical details: 
In appendixf^] we derive explicit expressions for the sym- 
metrized closed fcrmion loops of the FSM. In the follow- 
ing two appendices 151 and [Cl we carefully discuss the sym- 
metrized three-loop and the four-loop which are needed 
for the calculations in the main part of this work. Finally, 
in appendix |Pl we present a non-perturbative functional 
renormalization group flow equation for the irreducible 
polarization of the FSM. Although in this work we shall 
not attempt to further analyze this rather complicated 
integro-diffcrcntial equation, we use it in Sec. IIV Al to 
justify our self-consistency equation for n*(w, q). 



5 



II. RPA FOR THE FORWARD SCATTERING 
MODEL 

Because the RPA is exact for the TLM due to the 
vanishing of the symmetrized closed fermion loops with 
more than two external leg a 25 i 27 ' 32 i 33 , it seems at the first 
sight reasonable to use the RPA as a starting point of the 
perturbative calculation of the dynamic structure factor 
of the FSM. It turns out, however, that the RPA result 
for S(u), q) exhibits some unphysical features (see below) 
which arc related to the fact that the effect of interac- 
tions on the energy scale of the single-pair particle-hole 
continuum is not included in the RPA. A better starting 
point would be the so-called RPAE or "time-dependent 
Hartrec-Fock approximation" , because it takes the renor- 
malization of the single-pair particle-hole continuum ap- 
proximately into accoun t 13 i 15 . On the other hand, the 
RPA is sufficient to understand the relation between the 
mass-shell singularities and the expansion of the free po- 
larization in powers of 1/m, so that in this section, we 
shall carefully work out the spectral line-shape of the 
FSM using the simple RPA. The irreducible polarization 
is then approximated by the non-interacting one, 

n*(Q)«n (Q) = - / G (K)G (K + Q) , (2.1) 

JK 



where 



with 



G (AO = 



h 2 b 2 
2m 2m 



(2.2) 



(2.3) 



For f3 — ► oo and V — ► oo the integrations can be per- 
formed analytically, 

1 v e(-&) - e(-&+,) 



n„(Q) 



IU - £,k+q + £,k 



m 

— In 

irq 



tu + v F q+^ 



iu + v F q-^ 



(2.4) 



The corresponding RPA structure factor has been dis- 
cussed in Rcf. It consists of two contributions, 

Srpa(w,<7) = Z q 5(uj - uj q ) + Sg£ A (uj,q) , (2.5) 

where the first term represents the undamped ZS mode 
with weight 

,2 



(2.6) 



and energy 3 ^ 



v F \q\]l 1 + — coth 





q 




2k F 



v \q\ { 1 + 



5o(4 + 3 5o ) 



6a; 2 



2k F g 





FIG. 3: (Color online) Regime in the ui-q plane where 
5rpa(w, q) is finite. The shaded regime represents the single- 
pair particle-hole continuum, while the thick line corresponds 
to the ZS mode. For weak interaction jo C 1 the linear 
approximation u> q ~ Vo\q\ to the dispersion of the ZS mode 
crosses the upper boundary of the single-pair continuum at 
q ss gok F /2. However, the non-linear corrections to the ZS 
dispersion (|2.7|l are such that it never intersects the single-pair 
continuum, so that there is no Landau damping. 



The dimcnsionlcss function 



W 



kpgo 



sinh 



\ k F g J 



(2.8) 



can be identified with the relative contribution of the ZS 
peak to the /-sum rulei 2 - 



cLujuj S '(cj , q) = 



VFq 
2ir 



(2.9) 



The second part S^ A (ui,q) in Eq. (|2.5p represents the 
incoherent continuum due to excitations involving a sin- 
gle particle-hole pair (single-pair continuum)^. Because 
the ZS mode never touches the single-pair continuum, 
there is no Landau damping and within RPA the ZS 
mode is undamped. The damping of the ZS mode is 
due to excitations involving more than a single particle- 
hole pair (multi-pair excitations), which are neglected in 
RPA. The regime in the ui-q plane where S , rpa(o;, q) is 
finite is shown in Fig. [3J The corresponding qualitative 
shape of <Srpa(w, q) for fixed q -C k F is shown in Fig. [4] 
In the limit go — ■+ the ZS mode disappears and the inco- 
herent part S^p A (uj,q) reduces to the dynamic structure 
factor of the free Fermi gas, which in one dimension for 
q < 2k F is simply a box-function of width q 2 /m centered 
around i>_f|<?|) 



S (w,q) 



lim Srpa(w,<7) 

9o ->0 



2nq 



2m 



\U - v F \q\ 



(2.10) 
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(q/k^ 
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q/k F 



continuum 



1 



V /V F 



(0/v F q 



FIG. 4: (Color online) Schematic behavior of Srpa(w,?) for 
fixed q <C &f as a function of for q/Zcf <S 50 -C l.In this 
regime, the distance between the upper edge of the single- 
pair particle-hole continuum and the position of the ZS peak 
(indicated by a thick arrow) is much larger than the width of 
the particle-hole continuum. 



For finite go, the shape of 5ffp A (u>, q) is modified as shown 
quantitatively in Fig. 1 of Ref. [12| . The small shaded 
hump in Fig. [4] represents schematically the incoherent 
part of 5rpa(w, q) for finite go- For \q\/k F <C go, the 
relative weight of the single-pair continuum is negligi- 
bly small, so that the ZS peak carries most of the spec- 
tral weight. For example, the relative contribution of 
the single-pair continuum to the /-sum rule vanishes as 
(q/g k F ) 2 < I- 

It is instructive to see which features of SppA^q) 
are recovered if we expand the inverse non-interacting 
polarization IIq 1 (Q) in powers of the inverse mass wT 1 . 
Therefore we introduce the dimcnsionlcss variables, 



IU1 

W = — , P 
v F q 



q 

2k F 



and rewrite Eq. (|2.4[) as 

H.o(iu,q) = v fl (iy,p), 
with the dimcnsionlcss function 



(2.11) 



(2.12) 



rio(iy,p) = 7T ln 
2p 

4p 



iy + 1 + p 



iy + 1 — p 

'y 2 + {l+P? 



y 2 + (l-pf 



(2.13) 



For an interaction with momentum-transfer cutoff qo <C 
k F the relevant dimensionless momenta satisfy \p\ <C 1, 
so that we expand II^" 1 (z?/,p) in powers of p. From 
Eq. ([2T3]) we find 

n 2 1 - 3i/ 2 

uzHiv,p) = 1 + y 2 tt^4 + °(p 4 )- ( 2 - 14 ) 

o 1 + y* 

For later reference, we note that the correction of order 



p 2 in Eq. (|2.14[) can be written as 



P 2 1 - 3y 2 2 

= P 



3 (1 + y 2 ) 



1 - 



= P 



2p 2 



3(1 + y 2 ) 
1 



1 



1 - iy l + iy 



(2.15) 



For p — > we recover the result for linearized dispersion, 
limfi^iy.p) = tl^iw) = 1 + y 2 , (2.16) 

which yields the dynamic structure factor of the TLM 
given in Eq. (|1.2p . However, after analytic continua- 
tion to real frequencies iy — > x + iO = + iO, the 
correction term of order p 2 in the expansion (]2 . 14[) is 
singular on the mass-shell |cj| = vfM- Although in 
the non-interacting limit we know that this mass-shell 
singularity has been artificially generated by expanding 
the logarithm in Eq. (|2.13|) . it is not at all obvious how 
to regularize a similar singularity if it is encountered in 
the interacting system. This is the reason why a for- 
mal expansion in powers of the band curvature 1/ra us- 
ing either a purely fermionic approach^ or conventional 
bosonization^ is not reliable close to the mass-shell after 
analytic continuation. Fortunately, within the functional 
bosonization approach used in this work this problem 
does not arise, because the effective expansion parameter 
in functional bosonization is not 1/m, but the combina- 
tion goqo/(mv F ), see Refs. [25||27| . In particular, in the 
non-interacting limit functional bosonization yields the 
exact structure factor of the free Fermi gas with quadratic 
dispersion, containing all orders in 1/m. 

It is instructive to examine the RPA dynamic structure 
factor if we nevertheless use the expansion (|2.14[) for the 
non-interacting polarization. Then we obtain after ana- 
lytic continuation iy — > x -MO = u>/ {v F q) + iO for small 
kl < 9ok F , 



<Srpa(^,<7) ~ — Im 



1 



g + U^{x + iQ,p) 



= Z+5{u J -Cu+) + Z-8{La-u-) 1 (2.17) 

where and uj^ reduce for small q to the correspond- 
ing expressions Z q and Vo\q\ for linear dispersion [see 
Eq. (|1.4|) ]. and the weight and dispersion of the other 
mode Q~ is for \q\ -C k F go, 



z: 



v F \q\ 



2M 
3tt 



1 - 



-1 2 



2k F g 



3g \2k F 



(2.18) 



(2.19) 



This peak is associated with the incoherent part 
S^p A (uj,q) of the dynamic structure factor discussed 
above, which in the approximation (|2.14[) is replaced by a 
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single peak with the same weight. From Eqs. (|2.18[) and 
(|1.4p one easily verifies that for \q\ <C gok F the relative 
weight of the peak associated with the incoherent part is 
indeed small, 



4z 
3 



2k F g 



v F q 



foQo 



(2.20) 



where we have used po = qo/(2k F ), see Eqs. (|1.14[) . 
Hence, for \q\ <C gok F most of the weight of S , rpa(w,<z) 
is carried by the ZS mode uj q f=a Vo\q\, so that the inco- 
herent part corresponding to the mode ui~ ~ VF\q\ can 
be neglected^. Note that the limits q — ► and go — ► 
do not commute and that only for \q\/(2k F ) <C go the 
weight of the mode uj~ can be neglected. 

Mathematically, the second peak in Eq. (|2.17p is due to 
the pole arising from the term of order p 2 in the expansion 
()2.14p of the inverse free polarization. Although after an- 
alytic continuation iy — > x + iO this term is singular for 
x = 1 , we know from the exact result (|2.13|) how this sin- 
gularity should be regularized: we simply should smooth 
out the corresponding (^-function peak over an interval of 
width w g ex q 2 /m. In fact, we can self-consistently cal- 
culate w g by noting that after analytic continuation the 
singular term in the expansion (|2.14p gives rise to the 
following formally infinite imaginary part of the inverse 
non-interacting polarization, 



lmtl 1 (x+i0,p) = -T (x,p) 



2- 



P 2 [5(l-x)-6{l+x)}. 

(2.21) 

Ignoring the renormalization arising from the (singu- 
lar) real part of fl ~ 1 (x + i0,p) and approximating the 
resulting dynamic structure factor in this regime by a 
Lorentzian centered at u> = vfM, we find for the full 
width at half maximum in the limit go <C 1, 



VF\q\ 



T (l,p = q/(2k F )). 



(2.22) 



To obtain a self-consistent estimate for w q we follow 
Samokhin 1 and regularize the singularity in To(l,p) 
by replacing S(uj = 0) by the height of a normalized 
Lorentzian of width w g on resonance, 



S { x - l)U=a = v F \q\ S(uj - v F \q\ 



v F \q\ 



Hence, our self-consistent regularization is 

2p 2 v F \q\ 



3w 



(2.23) 



(2.24) 



Substituting this into Eq. (|2.22p wc obtain the self- 
consistency equation 



1 ( q Y (vFq) 2 



3 V 2k F ) w, 



(2.25) 



which leads to the following estimate for the width of the 
single-pair particle-hole continuum, 



1 



2-JZ' 



(2.26) 



Of course, it is now known^iii that the shape of 
the single-pair continuum cannot be approximated by a 
Lorentzian, but the order of magnitude of its width ob- 
tained within the above regularization is correct for suffi- 
ciently small q. Hence, the mass-shell singularity arising 
after analytic continuation iy — > x + iO in the expansion 
of the inverse non-interacting polarization (|2.14p in pow- 
ers of p — q/{2k F ) is simply related to the single-pair 
particle-hole continuum. This singularity can be regu- 
larized by smearing out the 5-function in the imaginary 
part over a finite interval of width w q cx q 2 jm. However, 
the width w q should not be confused with the damping 
of the ZS mode, which remains sharp within RPA. 

In order to obtain the ZS damping, one should cal- 
culate interaction corrections to the irreducible polar- 
ization. The finite overlap between the continuum due 
to particle-hole excitations involving more than a sin- 
gle particle-hole pair (multi-pair excitations) then de- 
termines the ZS damping. In three dimensions general 
phase space arguments^ 6 - imply that the resulting damp- 
ing is very small. In one dimension, an argument due to 
Teber— suggests that the damping of any acoustic collec- 
tive mode which overlaps with the two-pair continuum 
should vanish as q 3 for small q. However, for this argu- 
ment to be valid, one should self-consistently calculate 
the renormalized energy of the ZS mode and show that 
it is immersed in the multi-pair continuum. This has nei- 
ther been done in our previous worki^, nor in the work 
by Pustilnik et where the renormalization of the ZS 
velocity has been ignored. 

In this work we shall carefully examine all correc- 
tions to the RPA to second order in an expansion in 
powers of the small parameter po = qo/(2k F ), which 
is naturally generated using the functional bosonization 
approac h 25 ' 27 . Most importantly, our approach does not 
suffer from the mass-shell singularities discussed above. 
Moreover, we shall show that the distinction between 
the ZS energy vo\q\ and the energy scale v F \q\ associ- 
ated with the single-pair continuum shown schematically 
in Figs. [3] and [4] is an unphysical artefact of the RPA, 
which disappears once the corrections to the RPA are 
self-consistently taken into account. 



III. FUNCTIONAL BOSONIZATION 

In this section wc outline the functional bosoniza- 
tion approac h 25 ' 26 ' 27 which we then use in Sec. IIVI to 
calculate the dynamic structure factor. In contrast to 
Refs. 25.26,27], wc shall here keep track of Hartrce cor- 
rections to the fermionic self-energy, because these cor- 
rections contribute to the renormalization of the ZS ve- 
locity. 
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Decoupling the density-density interaction in Eq. (|1.8[) 
by means of a real Hubbard-Stratonovich field </>, the ra- 
tio of the partition functions with and without interaction 
can be written as 



fD[c, 



c,0e 



-So[E,c]-So[0]-Si[c,c,<? 



/ V[c, c, cj)}e- s °l 



c,c\ — Sa[ 



(3.1) 



where the free fermionic action So[c, c] is given in 
Eq. (|1.9[) . the free bosonic part is 



So[<f>] 



and the Fcrmi-Bosc interaction is 



Si [c,c, 



1 I / c k+qck(, 

IQ JK 



(3.2) 



(3.3) 



The fermionic part of the action in the numerator of 
Eq. (|3.ip can be written as 

S [c, c] +5i[c, c, cj)} =- / c k \G~ x \kk'Ck', 

J k Jk' 

(3.4) 

where the infinite matrix G -1 is defined by 

[G~ 1 \kk> = S K , K >[iuj -eje+ft]- i(f>K~K>- (3.5) 

At finite density the field </>a has a non-zero expectation 
value, 



= -*<$Q,O0 + 



(3.6) 



Here the <5-symbol is for finite /3 and V (where the com- 
ponents of Q = (u>,q) are discrete) given by Sq^q = 
PV5 a fi$q,o, which reduces to (2ir) 2 5(iD)5(q) for — > oo 
and V — ► oo. We fix the real constant </> from the re- 
quirement that the effective action ^efff^] of the 0-field, 
which is obtained by integrating over the fermionic fields 
in Eq. (|3.1[) . does not contain a term linear in the fluctu- 
ation 5<f>Q. To do this, we define the matrix G^ 1 which 
includes the self-energy correction due to the vacuum ex- 
pectation value <j), 



[G \kk' 



5 K .K'[ iuj - £ k - <j> + p], (3.7) 



and write 



with 



G 1 = Gn 



V, 



[V] A 'A'' = i$0K-K> 



(3. 



(3.9) 



Integrating in Eq. (|3.ip over the fermion fields we obtain 
the formally exact expression 



(3.10) 



where f2i — Oo is the change of the grand canonical poten- 
tial due to the vacuum expectation value ignoring fluc- 
tuations, 

fli - n = ^TrHGo^Go 1 ^ = 0)] - V^-. (3.11) 

The effective action for the fluctuations of the bosonic 
field is 



S e s[5(/>] = S [<j)Q -> -i5 Qfl (j> + 6<Pq] - /3V 
-Trln[l - G V] 

z Jq 



2/o 



Tr[G V] r 



(3.12) 



71=1 



We now fix the vaccum expectation value <j> from the 
saddle point condition 



or equivalently 



-Vf + V Po = 0, 
Jo 



4> = fopo- 



(3.13) 



(3.14) 



Here, po is the density and Gq(K) is the fermionic Green 
function in self-consistent Hartree approximation, where 



Po 



= f Go (JiQ = i V 0(a« - e fc - /opo) , (3-15) 



Go (K) 



1 



iu - e k - foPo + P 



(3.16) 



Note that Eq. (|3~T6)) agrees with Eq. (f2T2|) if we take into 
account that within self-consistent Hartree approxima- 
tion the Fermi momentum kp is defined via 



kp 
2m 



M - fopo- 



(3.17) 



Eq. (|3.13p guarantees that the terms linear in the fluctua- 
tions 6(f>Q in Eq. (|3.12|) cancel, so that our final result for 
the effective action for the fluctuations of the Hubbard- 
Stratonovich field is 



SesM = r / f- 



y\ Tr[G V] T 

n=2 



(3.18) 



= S 2 [6<i>]+Sto t [6<i>], 
with the Gaussian part given by 

W] = \l If- 1 + n o (Q)]<S0- Q ^ Q) (3.19) 
z Jo 




P(l) 




P(3) 



FIG. 5: Boson vertex with n external legs in the inter- 
action part Sint[<$0] of the bosonized effective action, see 
Eq. (|3,20p . The arrows denote the fermionic Green func- 
tions Go(K) within self-consistent Hartree approximation, see 
Eq. (|3.16[1 ■ The sum is taken over the n\ permutations of the 
labels of the external legs. For linearized energy dispersion 
all symmetrized closed fermion loops with more than two ex- 
ternal legs vanish. 




and the interaction part 

oc _^ 



Si, 



n=3 

r^(Qi, 



3 Qi + ...+Q„,0 



The vertices r£ ra) (Qi, 
lot 
2]), 

i Qn) 



symmetrized closed fermion loops L^(Q 
fined in Eq. 



,Qn)5<t>Qr---WQ n - (3-20) 

Q n ) are proportional to the 
■ , Qn) de- 



FIG. 6: Diagrams arising in the perturbative expansion of 
the irreducible polarization to second order in the RPA inter- 
action. The shaded circles represent the vertices of S c s[S^>], 
which are related to symmetrized closed fermion loops as 
defined in Fig. [S] Diagram (a) is equivalent to the three 
fermionic diagrams shown in Fig. [2] Diagram (b) is the so- 
called Aslamazov-Larkin diagram, while diagram (c) can be 
viewed as a higher order self-energy correction which renor- 
malizes the relation between density and chemical potential. 
Diagram (d) involving two bosonic loops and the symmetrized 
fermionic six- loop is of fourth order in po — qo/(2kp) and can 
be neglected to order p%. 



lt>(Qi 



*>-i)i4 n) (- 



A graphical representation of Iq' (Qi, 
in Fig. [5j In appendix [5] we give explicit expressions 
for the symmetrized n-loops of the FS M 28 i 29 i 30 and 
show that Lg oc (1/m)™ -2 to leading order in 1/m. 
Moreover, in appendices [B] and [C] we carefully discuss 
the properties of the loops Lg (Qi, Q2, — Qi — Q2) and 

(Qij Q27 — Qi, — Q2) for the special combinations of 
momenta needed in this work. 

The exact irreducible polarization can now be ob- 
tained from the fluctuation propagator of the Hubbard- 
Stratonovich field, 



Wick theorem. The RPA interaction thereby plays 
■Qi, . , . , — Qn). ro i e f thfi Gaussian propagator, so that we nat- 

(3.21) urally obtain an expansion in powers of the RPA in- 
teraction. As discussed in appendix \^ the vertices 
Q n ) is shown T (n) (Q 1; , , . ; Q n ) witn n external bosonic legs are propor- 



/D[^]e- s "«W] 
1 



J Q+Q'fl-—i 



fc L + ih(Q)' 



(3.22) 



where the effective action S e s [54>] is defined in Eq. (|3.18l) . 
Within the Gaussian approximation this reduces to the 
RPA interaction 



1 



5q+q ''W + MQ) 

£q+q',o/rpa(Q)- 



(3.23) 



The corrections to the RPA can now be calculated sys- 
tematically in powers of the interaction 5i n t using the 



tional to (1/m)™ -2 , so that the underlying small parame- 
ter of our expansion is the product of the RPA interaction 
and 1/m. In particular, in the limit of vanishing inter- 
action we recover the exact non-interacting polarization 
IIo(ia;,g) for quadratic energy dispersion. Our approach 
based on functional bosonization is therefore funda- 
mentally different from conventional bosonizatio n 1 ' 7 ' 14 , 
where the quadratic term in the energy dispersion gives 
rise to a cubic vertex proportional to 1/m which has to 
be re-summed to infinite order to recover the correct non- 
interacting polarization. 



IV. CALCULATION OF S(u), q) USING 
FUNCTIONAL BOSONIZATION 

A. One-loop self-consistency equation for II, (Q) 

It is now straightforward to expand the irreducible po- 
larization in powers of the RPA interaction, which is 
the Gaussian propagator of our boson field 5<pg. The 
diagrams contributing to II, (Q) up to second order in 
the RPA interaction are shown in Fig. [5] The relevant 
dimensionlcss parameter for this expansion is the ratio 
Po = 9o/(2&f) ; because by assumption the range of the 
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interaction in momentum space has a cutoff qo <C kp 
so that each additional bosonic loop integration gives 
rise to a factor of p\. Hence, the diagram (d) involv- 
ing two bosonic loops of is of order p^, because according 
to Eq. (|A17[) the symmetrized fermion loop with six ex- 
ternal legs is proportional to 1/m 4 and the two bosonic 



loop integrations generate a factor of q$. Because the 
other three diagrams are proportional to J5q, it is consis- 
tent to neglect diagram (d) as long as we retain all terms 
up to order Pq. Evaluating the diagrams (a)-(c) in Fig. [5] 
we obtain the following expression for the irreducible po- 
larization, 



n.(Q) « no(Q) - i^/RPA(go|64 4) (o^-Q^Q I -Q)+ 4 /RPA(o)4 3) (Q I -Q I o)4 3) (Q',-Q',o) 



+4/rpa(Q + Q')Lf (-Q, Q + Q', -Q')Lf (Q\ -Q - Q', Q) \. (4.1) 



r 



The properties of the symmetrized three- and four-loops 
appearing in this expression are discussed in detail in ap- 
pendices |B] and [Cj It turns out, however, that in order to 
cure the unphysical features of the RPA discussed at the 
end of Sec. HI] (in particular, within RPA the energy scale 
Vf\q\ of the single-pair continuum erroneously involves 
the bare Fermi velocity), we should self-consistently dress 
the Gaussian propagator /rpa(Q) in Eq. (|4.1|) by self- 
energy corrections. Formally, this amounts to replacing 
the RPA interaction by the exact effective interaction, 



/rpa(Q) f*(Q) 



i + /,n.(Q) 



(4.2) 



In appendix [D] we justify this procedure using a func- 
tional renormalization group approac h 36 ! 37 . With this 
substitution, Eq. (|4.1[) becomes a complicated integral 
equation for the irreducible polarization, which cannot 
be solved analytically. Fortunately, this integral equation 
can again be simplified by noting that on the right-hand 
side it is not necessary to retain the full Q-dependence 
of n*(Q), but to keep only those terms which contribute 
to the self-consistent renormalization of the ZS velocity. 
To explain this, let us introduce again the dimcnsionlcss 
variables iy = iuj/(vFq) and p = q/(2kp) and define the 
dimcnsionlcss irreducible polarization 



n*(zu;,g) = v a H.*(iy,p). 



(4.3) 



The corresponding dimcnsionlcss effective interaction is 
then 



9 P 



1 + g p IL*(iy,p) 



(4.4) 



where g p = vof q=p k F , see also Eqs. (|2.11l) and (I2.12[) . 
The dynamic structure factor can then be written as 



S(u),q) = — Im 

7T 



1 



Illl 



/ 9 + n7> + io,<z) 
i 

g p + Tl* l {x + iQ,p) 



(4.5) 



where x = cu/(vpq). For our purpose it is now sufficient 
to approximate the dimensionless inverse irreducible po- 
larization by 



ti- 1 {iy 1 p) = Z l + Z 2 y 2 



(4.6) 



where the dimcnsionlcss renormalization factors Z\ and 
Zi should be determined as a function of the interac- 
tion such that the approximation (|4.6[) yields the true ZS 
velocity v. Within RPA, where the non-linear terms in 
the energy dispersion do not renormalize the ZS veloc- 
ity, the irreducible polarization is approximated by the 
non-interacting one, so that Z\ = Zi = 1. If we approx- 
imate the inverse polarization in Eq. (|4. 5[) by Eq. (|4. 6|) 
we obtain for lo > and q — > 0, 



VF\q\ 



8{uj - v\q\), 



2izvZ2 

where the renormalized ZS velocity is 



v 

v F 



Zn 



X 



with renormalized coupling constant 



go + Z 1 



- 1. 



(4.7) 



(4.8) 



(4.9) 



In order to avoid the unphysical splitting of the spectral 
weight in S((J, q) (as discussed at the end of Sec. lU this 
is an artefact of the RPA) it is crucial that the true ZS 
velocity v appears in the bosonic propagators. There- 
fore, a naive expansion in powers of the RPA interaction 
is not sufficient. However, we may further reduce the 
complexity of the calculation by noting that Eq. (|4.7p 
still contains the correct velocity if we set Z-2 — * 1 in 
the prefactor. Within this approximation, the velocity 
renormalization implied by Eq. (|4.6[) can be simply taken 
into account via a re-definition of the coupling constant, 
go — ► g. It is therefore sufficient to replace the RPA in- 
teraction in Eq. (|4.1[) by an effective interaction of the 
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same form but with a rcnormalizcd effective coupling g 
instead of go, which should be chosen such that all inter- 
action corrections to the ZS velocity are self-consistently 
taken into account. Note that Schonhammer— has re- 
cently shown that within the so-called RPAE approxima- 
tion (which amounts to solving the Bethe-Salpeter equa- 
tion with the bare interaction as irreducible vertex) the 
relative position of the collective mode energy and the en- 
ergy of the single-pair particle-hole continuum is different 
from the RPA prediction for the FSM: In RPAE the ZS 
mode lies above the non-interacting single particle-hole 
continuum which (erroneously) appears in RPA, but be- 
low the Hartree-Fock particle-hole continuum. This sug- 
gests that in order to obtain a correct estimate of the ZS 
damping, it is necessary to calculate the location of the 
ZS energy self-consistently. 

In field-theoretical language the constants Z\ and Zi 
are counter-terms which guarantee that our Gaussian 
propagator depends on the true ZS velocity. In Sec. [V] 
we shall explicitly calculate the factors Z\, Zi and the 
corresponding rcnormalizcd ZS velocity v to second order 
in our small parameter po- A similar procedure is neces- 
sary to self-consistently calculate the true Fermi surface 
of an interacting Fermi syste m 38 ' 39 . The expansion of the 
modified dimensionless interaction g p for small p is then 

~9v = 9+\gy +0(p 4 ), (4.10) 



where 

g' ' = (2fc F )Wo = sign/£ — . (4.11) 

np c 

In this approximation, our dimensionless effective inter- 
action is 

f*(iy,P)*f g (w,P) = z r , (4-12) 
1 + g p U {iy,p) 

which differs from the RPA interaction, because the func- 
tion g p includes the renormalization of the ZS velocity 
due to fluctuations beyond the RPA. 

Collecting all terms, our final result for the dimension- 
less irreducible polarization to one bosonic loop can be 
written as 

&»(<!/, p) « tl (iy,p) + Bitty, p) + ti 2 {iy,p), (4.13) 

where the non-interacting polarization is given in 
Eq. (|2.13p . and the subscripts indicate the powers of g p . 
The term ITi (iy,p) corresponding to diagram (a) in Fig. [6] 
can be written as 



/OO pOO J I 

dp'\ P '\ / ^J 9 (iy',p')if(iy,p,iy',p'), (4.14) 
- oo J — oo 

where the dimensionless symmetrized four- loop L^\iy,p,iy' ,p') is defined in Eq. (|C9|) . The term fl2(iy,p) involving 
two powers of the effective interaction is of the form 

U 2 (iy,p) = IL^(iy,p) + tLf(iy,p), (4.15) 
where the contribution from the Aslamasov-Larkin (AL) diagram in Fig. [5] (b) is 

ftf^p) = - I" dp'\p'\ ^ ^f g (iy',p')f g ( iyP + iy ', P ' ,p + p') [4 3 W, i2 /,;/)] 2 , (4.16) 

J-oo J-oo 27T \ P + P J 

and the contribution from the Hartree diagram in Fig. [5] (c) can be written as 

/oo />oo 7 / 

d P '\p'\ / ^~fSy',p')L { s\iy'iP\iy'-p')- (4.17) 
-oo J-oo ^ 7r 

I 



Here, the dimensionless symmetrized three-loop 
L^g\iy,p,iy' ,p') is defined in Eq. (|B3jl . The parameters 
Z\ and Zi hidden in the effective interaction f g (iy,p) 
should be determined self-consistently by evaluating 
Eqs. (|4.13H4.17p and demanding that the resulting 
rcnormalizcd ZS velocity is consistent with the result 
obtained from Eq. (|4.6| . We emphasize again that 



the above expression for II, (Q) is not based on an 
expansion in powers of 1/m: all functions appearing in 
Eqs. ((4~T4H4~T7|) depend on 1/ m in a rather complicated 
non-linear manner. 
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B. Approximation A: neglecting 1/m-corrections to 
Ho(Q) in loop integrations 

Eqs. (|4.14I - |4~T7|) are still too complicated to admit an 
analytic evaluation. In order to explicitly calculate the 
dynamic structure factor without resorting to elaborate 
numerics, wc shall further simplify the above expressions 
by making the following approximation A: We replace 
the non-interacting polarization Hq(Q) appearing in the 
effective interaction and the symmetrized closed fermion 
loops on the right-hand sides of Eqs. |^.i^|[4TT7} ) by its 
asymptotic limit for small momenta given in Eq. \2.16\ ). 
Keeping in mind that in one dimension the closed fermion 
loops with n > 2 external legs can all be expressed in 
terms of Tlo(Q), the symmetrized three- and four-loops 
are then approximated by Eqs. (|B8[) and (|C14[) . For con- 
sistency, we should also expand the dimensionless free 
polarization Ho(iy,p) on the right-hand side of Eq. (|4. 13|) 
to second order in p, see Eq. (|2.14p . Wc shall argue below 
that the above approximation A is not sufficient to cal- 
culate the line-shape of the dynamic structure factor for 
momenta q < q c = l/(m|/o |) [see Eq. (|1.16[) ]. because in 
this regime the spectral line-shape is dominated by the 
terms neglected in approximation A. On the other hand, 
for q > q c the line-shape of S(uj 7 q) is essentially deter- 
mined by the quadratic term in the expansion of f q for 
small q, so that in this regime A is justified. 

It turns out that with this simplification the y'- 
integrations in Eqs. (|4.14[) . (|4. 16|) and (|4.17[) can be done 
analytically for general g p using the method of residues. 



The form of Eq. (|4. 5|) suggests that it is natural to expand 
the inverse irreducible polarization in powers of p and p Q . 
This procedure can be formally justified within functional 
bosonizatio n 12 i 27 , where the interaction corrections to 
the inverse irreducible polarization play the role of the 
sclf-encrgy corrections in the effective bosonized theory. 
But it is usually better to expand the self-energy rather 
than the Green function in powers of the relevant small 
parameter, because the direct expansion of the Green 
function usually leads to unphysical singularities. Using 
Eqs. (|2.14|) and (|4. 14114.17) we obtain for the expansion 
of the inverse irreducible polarization to order p 2 ,, 

n* {iy,p) = 1 + 2/ --5- 1,2 

3 1 + y £ 

- (l + y 2 ) 2 n 2 (iy,p) + O(p 3 ), (4.18) 

where p is assumed to be smaller than the dimensionless 
momentum-transfer cutoff po = qo/(2kp). It is conve- 
nient to introduce the notation 



x P = v 71 +S P , (4.19a) 
a p = Xp + 1 = y/l + g p + 1, (4.19b) 
bp = x p -l = y/1 + g p - 1, (4.19c) 

so that a p b p = g p . The contribution involving the sym- 
metrized four-loop can then be written as 



with 



(l + yzyn^p) = Re / dp' 



\p'\p' A F 1 (iy 1 p')+p' 2 p 2 F 2 (iy,p')+p i F 3 (iy,p') 



Kp' 2 



(1 + iy) 2 p 2 ][b 2 ,p' 2 - (1 - iy) 2 p 2 



-p' 2 g pl {l+iyf 



p(l - iy) 



\p+p'\ 



x 2 ,p' 2 



[(1 - iy)p + p 



'12 



+ (p' - -p') k(4.20) 



Fi{iy,p) 
F2(iy,p) 
F 3 (iy,p) 



4ff p (a; p + iy) 2 + g 2 



8x p 
1-iy 



9 P 



4x p -g p -(2 + x p -^ )(l+y 2 ) 



-(1 + iyf + ~g p (2 -y 2 + y 4 ) - 4b p iy(l - y 2 )] - 2b 2 p x p \±^(fi - 6y 2 - y 4 ) 

1 zy 



4^(1 + y 2 ) 



l + y 2 (l + y 2 ) 2 



(4.21) 
(4.22) 
(4.23) 



Both functions Fi(iy,p) and F2(iy,p) contain a singu- 
lar term proportional to (1 — iy)" 1 , which after analytic 
continuation give rise to a mass-shell singularity at the 
energies ±VFq associated with the bare Fermi velocity. 
Fortunately, these singularities cancel when Eq. (|4.20[) is 
combined with the corresponding contributions from the 



expansion of n^" 1 (iy,p) in Eq. (|4. 1 8f> and from the AL 
diagram given in Eq. (|4.30[) below. To show this explic- 
itly, it is useful to isolate the singular term in Eqs. (|4.21|) 
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and (|4.22p by setting 



F\(iy,p) and F 2 (iy,p) arc now analytic functions of y, 



Fi{iy,p) = -+F 1 (iy,p), 



1-iy 

F 2 (iy,p) = -uv p u, p 



(4.24) 



*olx J\ + iy) +F 2 {iy,p). (4.25) 



1 — iy 



Fi(iy,p) = ^9p(x P + iy) 2 - ~g 2 p 



4x p + g p + (2 + x p - ^)(1 + y 2 ) 



F 2 (iy,p) = 9 P -{l + iyf + g P (2-y 2 + y A ) + Ab p iy{\ + 2iy + y 2 ) 



+ 2b 2 p {l + iy) [x p (l + iy)(l + y 2 ) - 4iy] . 



(4.26) 



(4.27) 



Eq. (|4.20p can then be written as 

poo 

-(l + 2/ 2 ) 2 ni(«/,p) = Re / dp 



2xafW- n t> l 0O , l \p'\p' 4 F 1 ( l y,p')+p' 2 p 2 F 2 ( t y,p')+p 4 F 3 ( t y,p') i 8\p'\ 8\p'\p 2 (l - iy) 



V Wp'V a - (1 + iy) 2 P 2 ][b 2 p ,p 12 - (1 - iy) 2 p 2 } 1 - iy b 2 ,p' 2 - (1 - iy) 2 p 2 



p' 2 ~g p ,{l + iyY 



2p> 



\P'+P\ 



p(l - iy) J x 2 ,p' 2 - [p' + (1 - iy)p) 



+ (p'^-p')\- (4.28) 



Next, consider the contribution TI 2 h (iy,p) from the Aslamasov-Larkin diagram in Eq. (|4.16p . Adopting again 
approximation A, the symmetrized three-loop L s (iy,p,iy' ,p') is replaced by its limit L s {iy,iy' ,p/p') for 1/m — > 
given in Eq. (|B8[) . Then we obtain 



-(l + y 2 ) 2 fl^(iy,p) 



dp'\p'\g p > 

9p'+p 



dy' 



i-yy 1 -{y + y'Y-^ 



2tt 



l + y>- 



x 2 p , + y' 2 



( py+p'y' \ ' 

V p+p' J 



2 , I py+p'y' 
^p'+p ^ \ p+p' 



2n ' 



(4.29) 



The j/'-integration can now be carried out using the theorem of residues. The result can be cast into the following 
form, 



(l+y 2 ) 2 fT^ L (*y,p) = Re / dp'p 



i ,\p'+p\\ Sp'+p\p' +P\[(p' +p)(l + 2iyx p , +x 2 ,) -p(i i 2 + x 2 pl ) 



, Xj/ [{p 1 + p) 2 ~(x p >p' + iyp) ] [x 2 pl+p {p' + p) 2 - (x p ,p' + iyp) ] 



g P 'P 



p'(l + 2iyx pl+p + x 2 pl+p ) +p(y 2 + x 2 pl+p ) 



Xp'+p[p' 2 -{x P '+p{p' +p) - iyp) 2 ] [x 2 p >P 12 ~ {x P +p'{p + p') - iypf] 



2(p'+p) 



gp'+pW +p\{i + iyf 



p(l - iy) J x 2 pl+p {p' +p) 2 - \p'+p- (1 - iy)pf 



2p> 



gp'p'il + iyf 



p(l - iy) J x 2 ,p' 2 - \p' + (1 - iy)p] 2 



+ (p->-p). (4.30) 



Finally, the contribution (|4. 1 T[) of the Hartree-type-of with 
diagram (c) in Fig.[6]is^ 



(l + y 2 ) 2 Tl 2 i (iy 1 p) = I H (l~y 2 ), (4.31) 



Ih = ~ 



29 



1 + .9 J, 

g 



i + g Jo 



dpp 



dpp 



Up 



(z p -l) 2 



- 1 



(4.32) 
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For 0-function cutoff this reduces to 



Ih = - 



Po9 
1 + 3 

while for Lorcntzian cutoff, 



1 + 2 

=1=-1 



Ih = 



via 

1 + 3 



9 



where we have used Eqs. (|2.14[) and (|2. 15[) to clearly ex- 
hibit the mass-shell singularity generated by the cxpan- 
(4 33) S10n °f the inverse free polarization. The dimensionless 
integral I{iy,p) can be written as 



(4.34) 



Combining all terms we obtain the following expansion 
of the inverse irreducible polarization to second order in 



Po> 



f£ l (iy,p) = l + y 2 + P 2 -%- 



1 — iy 1 + iy 



i(iy,p) 



dp'p' [J(iy,p,p') + J(-iy,p,p')] 



+ijy(l — y ) + I(iy,p) + 0(pq), (4.35) where the complex function J(iy,p,p') is given by 



(4.36) 



J(iy,p,p') 



p' 4 F 1 (iy,pi)+p> 2 p 2 F 2 (iy,p')+p 4 F 3 (iy 1 p') 8 



\l-iy) 



V Wp'P' 2 - (1 + iy) 2 P 2 } [bp/P 12 - (1 - iy) 2 P 2 } l — iy b 2 ,p' 2 - (1 - iy) 2 p 2 

r i 2 

|p'+p| f V+pIp'+pI (p / +p)(l + 2i3/ar p / + x 2 pl ) - p(y 2 + x 2 pl ) 



x p > [{p' + p) 2 - {Xp>p' + iyp) 2 ] [x 2 p , +p {p' + p) 2 - (VP' + iyp) 

2 



p'(l + 2iyx p , +p + x 2 pl+p )+p{y 2 + x 2 pl+p ) 



Xp>+p[p' 2 -(xp, + p{p' +p) - iyp) ] [x 2 ,p' 2 - (x p+ p,{p + p') - iyp) ] 



2{p'+p) 



9p'+p\p' ' +p\{l+iyf 



p(l - iy) J x 2 , +p {p' + p) 2 - \p> +p- (1 - iy)p] 



2p' 



1 



g p >p'(l + iyf 



p(l - iy) 1 x 2 ,p' 2 - \p' + (1 - iy)p] 2 
I 



+ (p-> -p). 



(4.37) 



Although it is not obvious from Eq. (|4.37[) . the function To proof this, it is sufficient to calculate the residues 
J(iy,p,p') vanishes as g 2 , for p' 3> po, so that the inte- 
gral ([4T36]) is ultraviolet convergent as long as g p vanishes R±(p) = 1(1 + x )I(x,p)] 
faster than l/p for p — > oo. 



C. Cancellation of the mass-shell singularities at 

lo — ±w_pg 

We now show that the mass-shell singularities at iy — > 
x = ±1 (corresponding to frequencies w = ±vf<i) arising 
from the expansion of the non-interacting polarization in 
Eq. (|4.35p are exactly cancelled by corresponding singu- 
larities in I(x,p), because for x — > ±1 the integral I(x,p) 
diverges as 



I{x,p) 



2p 2 1 
3 Itx 



x -> ±1. 



(4.38) 



i r°° 

- / dp'p' \im[(lTx)J(±x,p,p')}. (4.39) 



Using - 1 = g p we find from Eq. (|4.37|) . 

lim [(1 T ,)J( ± W)] = S- ^+d-fr'-*' 



80(b|-p')(l-p7N)- (4-40) 



Hence, 



i2±(p)=4 J P dp'p'(l-p'/\p\) = ^.. 



(4.41) 



which proofs Eq. (|4. 38)) . We conclude that the expansion 
(|4.35p of the inverse irreducible polarization to second 
order in p\ does not exhibit any mass-shell singularities 
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at frequencies u = ivpq corresponding to the excitation 
energy of non-interacting particle-hole pairs. This can- 
cellation also corrects the unphysical feature of the RPA 
that the single particle-hole pair continuum is centered 
at the energy Vp\q\ involving the bare Fermi velocity vf- 
It is convenient to explicitly cancel the mass-shell sin- 
gularities arising from the expansion of the free polariza- 
tion in Eq. (|4.35|) against the corresponding singularities 
in I(iy,p). Therefore we use the identity 



IT 
1 



1 



1 



1 — iy 1 + iy 

oo 

d P 'p'[J {iy lPl p') + J (-iy,P,P% (4.42) 



where 



1-iy 



p'+P+\p'+p\ 



2p 



+ (p ~P) 



(4.43) 



to write Eq. (|4.35|) as follows, 

fif^iw.p) = l + y 2 +P 2 + I H (l-y 2 ) + I(iy,p) 

+0(pg). (4.44) 

The integral I(iy,p) can again be written as 



i{iy,p) = 2 J d p'p' 



J(iy,p,p') + J(-w,P:p') 



with 



(4.45) 



J(iy,p,p') = J{iy,p,p') - Miy>PrP')- (4-46) 

We may now explicitly cancel the mass-shell singularities 
in the regularized integrand J(iy,p,p') and obtain after 
some algebra, 



J(iy,p,p') 



\i-iy) 



p> A F 1 (iy,p')+p' 2 p 2 F 2 (iy,p') + p ^F 3 (iy,p') | 

Xp> [a 2 p ,p' 2 - (1 + iy) 2 p 2 ] [b 2 p ,p' 2 - (1 - iy) 2 P 2 ] b 2 ,p' 2 - (1 - iy) 2 p 2 

9p'+p(p'+p) 2 \(P'+ P)(1 + liyXp, +x 2 )- p(y 2 + x 2 ) 



2x p > [(p' +p) 2 - (xp>p' + iyp) 2 } [x 2 p , +p {p' +p) 2 - {xp>p' + iyp) 
9 P '\P' +P\p' + 2iyxp> +p + x 2 pl+p ) +p(y 2 + x 2 pl+p ) 



2xp> +p [p n -{xp> + p(p' +p) - iyp) ] [x 2 pl p 12 - (x p+p '(p + p r ) - iyp) ] 

(p> + p) Wp' + p) - 4(1 - iy)p + ~g pl+p (p' + p) [*±B. ( 3 + iy) + 1(1 + iy) 2 



W+P\ 



X l'+ P (p' +p) 2 -\p' +P~ (!-*y)p] 2 

4(1 - iy)p + ~g p ,p' [£ (3 + iy) - |(1 + iy) 2 ] 



x 2 p ,p' 2 - \p' + (1 - iy)pf 



+ (p^ -p)- 



(4.47) 



V. INTERACTION WITH SHARP 
MOMENTUM-TRANSFER CUTOFF 

A. Explicit evaluation of the irreducible 
polarization 

In this section we assume that the dimensionlcss inter- 
action g p is of the form 

9p = 9oQ(po - \p\)- (5.1) 

In this case the p'-integration in Eq. (|4.36p is elementary 
and can be carried out exactly. Note that all derivatives 
of the interaction (|5.1[) vanish at p = so that /q = 0, 
which is certainly an unphysical feature of the O-function 
cutoff. The length q c defined in Eq. (|1.16[) is then for- 



mally infinite, so that the regime (|1.17[) does not exist. 
Although for such an interaction approximation A dis- 
cussed in Sec. lIVBl fi.e.. replacing Ilo(iy,p) rts flo(iy, 0) = 
[1 + y 2 ] -1 in loop integrations) is never justified, it is still 
instructive to evaluate Eq. (|4.44|) . because it allows us 
to explicitly see the partial cancellation between contri- 
butions arising from the first-order diagram in Fig. [5] (a) 
and the AL diagram in Fig. [6] (b). To clearly exhibit 
this cancellation, it is instructive to evaluate the contri- 
butions Hi(iy,p) (first-order in the effective interaction) 
and Il2(iy,p) (second order in the effective interaction) 
separately. Therefore, we specify g p = gO(po — |p|) in 
Eqs. (|4.36I4.37| and perform the p'-integration exactly. 
Recall that the effective coupling constant g is defined 
as a function of the bare coupling g via Eq. (|4.9p . The 
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p — > limits of the coefficients x p , a p and b p given in 
Eqs. (|4.19all4".19c[l are now denoted by 



6 



V 71 + .9, 
a?o + 1, 
a: - 1. 



(5.2a) 
(5.2b) 
(5.2c) 



After some tedious algebra we find that the contribution 
from the diagram (a) in Fig. [6] to the expansion (|4.18p 
can be written as 



Note that for small g, 



a = 2 + f - ^ + |! + 0( ff 4 ), (5.3) 
b = o-2=| - Y + f^ + °(3 4 )- ( 5 - 4 ) 



-(n-irrni(iy,p) = -pg 



2 6 2 (3 + x ) 



(2 + ff-A) 



2 J 2 1 - 3y 2 (2 + g) , , 4A 
P < t;^—^- + — — (4 ~9)-y 



3 1 + y 2 
J-A 



4 + -^t 



Re 



— (1 - zy)(x - «y) In ' 



+ 13 ( X ~ «y)(zo + In 



p 2 (l + iy)(x Q - iy) 
1 + iy 



(5.5) 



where we have defined 

A = 1 +g + y 2 = xl +y 2 . (5.6) 

If we neglect at this point the contribution Tl2(iy,p) in- 
volving two powers of the effective interaction, we re- 
cover from the imaginary part of Eq. (|5.5p our previous 
estimate^ 2 - for the damping of the ZS mode for q — > 

ir g 3 \q\ 3 

lq ~ T 1 2- ^ 5 - 7 ) 

8 Xoa 4 vpm z 

In view of the discussion at the end of Sec. [TT| this result 
should not be surprising: within our approximation the 



ZS mode is located at higher energy than the single-pair 
continuum and is immersed in the multi-pair continuum, 
whose spectral weight is generated by the logarithmic 
terms in Eq. (|5.5|) . The overlap of the multi-pair con- 
tinuum with the ZS mode leads to the q 3 -damping, in 
agreement with the arguments by Teber— . 

Unfortunately, the term in Eq. (|5.5[) which is respon- 
sible for the result (|5.7|) is exactly cancelled by a simi- 
lar term in — (1 + y 2 )!!^ (iy , p) . Explicitly carrying out 
the p'-integration in Eq. (|4.30[) and adding the contribu- 
tion (|4.3ip from the Hartree-type-of term, we obtain for 
\p\ < Po, 
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(l + y 2 ) A Tl 2 (iy 7 p) 



p2 °24 9{2 + g ~ A) 



+Po(po - A, [5(2 + 9) - 6(1 + |) - Ax 2 



(4 + g) 2 + 83(2 + g — A) , 5 2 A 



1 - 


3y 2 


3 (2 + 9) 


4 - 

4- 9+ — 
x - 


2Ar 


3(1 


f y 2 ) 


2x 2<? 





3.9 2 



12x n A 



16x^ 



In 



4 + .g-^- + -^+x (4- ff ) 
4 4xo 

4po(Po- b|)xg+p 2 A 



•'•() 



-Rc 



r(l - iy){xo - «2/)ln 



p 2 A 

Po(Po ~ \p\)a 2 + P 2 (l + w)(xo - iy) 



p 2 (l + iy)(x - iy) 



--jj( 1 -iy){xo + iy) In 



1 



iy 



xo + iy 



(5. 



Adding Eqs. (|5.5|) and Q5.8|) and rearranging terms, we obtain for the expansion (|4.18p of the inverse irreducible 
polarization for sharp momentum-transfer cutoff 

n* 1 ^,^) = 1 +Po9i + (1 +Po92)y 2 +Pa\p\[93 + 94y 2 ] 



b_ 

gxo 

(4 + 3 5 ) 2 , 5 2 A 



%- \ p- - 2 + — [8 + Ag - g 2 } + ^ [166 - Aga + g 2 (l + 3/x ) 
2 I 3xn QXo g 



6xoA 

<i + y 2 ) 



$4 

b 2 

2 ^-?-2Re 

a 6 Xo 



In 



4po(Po~ \p\)x 2 +p 2 A 



p 2 A 

(1 - iy)(x - iy) In 
I 



p a - \p\(xp - iy) 
p a + \p\(l + iy) 



(5.9) 



where 



b 2 




Q , 9x0 + 3" 






2xq 


L + 2^ + l_ 






3 2 
8 a 


5 3 
8 


-0(3 4 


)• 


(5.10a) 


b 2 

9 T 3 
zx 




V- 

8 y 


1 3 
4 5 + 


0(<? 4 ), 


(5.10b) 


6 2 


2^0 + ^0 


1 

= 8 g 


2 - -3 3 - 
32 y 


hO(g 4 ), (5.10c) 


b 2 

axo 




V- 

8 y 


5 3 
32 


-0(.9 4 )- 


(5.10d) 



Eq. (|5.9p has three important properties: 

• The logarithmic term in Eq. (|5.5p which is respon- 
sible for the g 3 -dependcncc of j q in Eq. (|5.7p . is 
exactly cancelled by a similar term with opposite 
sign arising from the AL diagram. 

• The mass-shell singularity at u> = ±vf<7 associated 
with the expansion of the free polarization n (w, q) 
in Eq. (|4.35|) has disappeared in Eq. (|5.5p , in agree- 
ment with our general considerations in Sec. IIV Cl 



• Eq. (|5.9p contains a term proportional to 1/A, 
which after analytic continuation gives rise to a 
mass-shell singularity at the physical energy u = 
±vq of the ZS mode. 

The mass-shell singularity at u> — ±vq is an artefact 
of the sharp momentum-transfer cutoff used in this sec- 
tion in combination with approximation A discussed in 
Sec. IIVBI In fact, we shall show in Sec. I VII that a more 
realistic interaction f q with finite /q does not lead to any 
mass-shell singularities, even if we still use approximation 
A to evaluate Eqs. (|4T4H4~T7)) . 



B. Renormalized ZS velocity 

To calculate the renormalized ZS velocity it is sufficient 
to set p — in Eq. (|5.9p . so that the problems related 
to the mass-shell singularity do not arise. Comparing 
Eq. (|5.9p at p = with the defining equation (|4.6p of the 
renormalization constants Z\ and Z2, we find to order 
Pi 

Zi = 1+pI Si , i = 1, 2, (5.11) 
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which are non-linear self-consistency equations for Z\ and 
Z 2 , because g\ and g 2 are defined in terms of the renor- 
malized coupling g = (go + Z\ — Z 2 )/Z 2l see Eq. (|4.9|) . 
However, keeping in mind that the difference g — go is 
proportional to p\ and that Eq. (|5.11[) is only valid to or- 
der Pq, we may ignore the self-consistency condition and 
set Z\ = Z 2 = 1 in the expressions for g\ and g 2 on the 
right-hand side of Eq. (|5.1ip . From Eq. (|4.8[) we then 
obtain for the rcnormalized ZS velocity, 



where vq == vp^l + go is the RPA result for the ZS 
velocity. A graph of the relative change of the ZS velocity 
as a function of the interaction strength g is shown in 
Fig. \7\ Obviously, even for large g and p\ = Oil) the 
correction to the RPA result «o never exceeds more than 
a few percent. 



v 

1>F 



where 



Zi + go 



Z 2 



(5.12) 



Ad hoc regularization of the mass-shell 
singularity and spectral line-shape 



.9 =9o -P 95, 



with 



.95 = x Q g 2 - gi = — 

Xr\ 



1 „2 ^ 3 



0(5 4 )- 



(5.13) 



(5.14) 



To order p\ we thus obtain for the energy of the ZS mode 



u q w v\q\, 



(5.15) 



Although for sharp momentum-transfer cutoff the dy- 
namic structure factor exhibits (within approximation A 
discussed in Sec. IIVB[) a mass-shell singularity at the 
ZS energy v\q\, it is nevertheless instructive to follow 
Samokhini and regularize the singularity by hand using 
the procedure outlined in Sec. [Til Because the natural 
scale for the momentum dependence is not 2fci? but the 
scale <7o set by the momentum-transfer cutoff, it is conve- 
nient to express the momentum dependence via q = q/qo- 
Setting p = poq and writing 



with renormalized ZS velocity, 



v F\/l + ga-plg5 



PQ7TJ2+O{po) 
Zi 



S(w, q) = — Im 

7T 



.go + n~ 1 (x + iO,g) 



(5.17) 



(5.16) we obtain on the imaginary frequency axis 



go+n^{iy,q) = A[l+p 2 {g 2 +g 4 \q\)] -p 2 g 6 \q\ 



Pal 



.97 + .98 mil 



4sg(l-|g|) 
<f A 



+ (.9 - A) 



a 3 xo 



Rc 



(1 - iy)(xo - iy) ln( — 
V a 



- \q\(x - iy) 



(5.18) 



where 



96 

.97 

98 

h 
hi 



x 2 g4 



93 



ax 



o 



b ( q 
g V 4 



9 



1 



16x1 



:9 



= -1 



16" 32 

ixo -Zgxo 
(l + 3x 2 ) 2 (4 + 3 5 ) 2 



9 3 + 0(g% 
l + 4g-g 2 ] = l 



TE 9 



1 3 8 (a b\ 1 2 11 o 4 . 



— q q 

6 y 12 y 



12a;o 



12x 



0(g% 



+ ^9+- A 9 2 + 0{g 3 ). 



0(g 3 



(5.19a) 
(5.19b) 
(5.19c) 
(5.19d) 
(5.19e) 
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FIG. 7: (Color online) Relative renormalization Av/(dqPo) = 
(v — vo)/(voPo) = — gs/(2xo) of the ZS velocity in units of Pq 
as a function of the interaction strength g, see Eq. (|5.16l l. 



From Eq. (|5.18p it is obvious that our functional 
bosonization approach yields a systematic expansion of 
the inverse irreducible polarization in powers of the small 
parameter po = qo/(2kp). Note that only ho and hi 
have finite limits for g — > 0, whereas the other couplings 
gi, . . . ,gs vanish at least as g 2 (the coupling _gg vanishes 
even as g 3 ). 

In the limit g — > Eq. (|5.18[) correctly reduces to 
the expansion of the non-interacting inverse polarization 
given in Eq. (|2.14p . However, the term hi /A generates a 
mass-shell singularity at the true collective mode energy 
w = ±vq. Fortunately, this singularity can be avoided if 
we use a more physical interaction whose Fourier trans- 
form fq is analytic for small q, as will be shown ex- 
plicitly in Sec. IVII In this subsection we shall sim- 
ply regularize the mass-shell singularity by hand using 
the self-consistent regularization procedure proposed by 
Samokhini, which we have already described in detail 
in Sec. [TTJ Repeating the steps leading from Eq. (|2.21|) 
to Eq. (|2.26p . we obtain from the self-consistent regu- 
larization of the singular term proportional to h\j A in 
Eq. (|5.18p the following estimate for the width of the ZS 
mode, 



sfhl(? 
2xq 2m 



(5.20) 



where we have factored out the corresponding estimate 
in the absence of interactions given in Eq. (|2.26|) , and the 
dimcnsionless factor Z w is given by 



I3hi 
4x1 
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is 



[l + 5] 3/4 ' 



(5.21) 



Note that Z w ~ 1 + ^g 2 + 0(g 3 ) for g -> 0, and 
Z w ~ jg 1 / 4 for g — > oo. A graph of Z w as a function 
of the interaction strength g is shown in Fig. [5] The es- 




FIG. 8: (Color online) Graph of the factor Z w defined 
in Eq. (|5,21[) , which estimates the interaction-induced rela- 
tive change of the width of ZS resonance for q <JC q c , see 
Eqs. (|5.20I5.21|) . 



timate (|5.20[) for the width of the ZS resonance on the 
frequency axis scales as q 2 , which is for small q much 
larger than our previous estimate (|5.7p based on the eval- 
uation of only the first-order diagram (a) in Fig. [5] The 
g 2 -scaling of the width of the ZS resonance has already 
been found by Samokhini and has been confirmed later 
in Refs.@lllg. However, the derivation of Eq. (j5T2TJ)) is 
based on a rather ad hoc regularization prescription of 
the mass-shell singularity in Eq. (|5.18p . which ignores in 
particular the divergent real part of the term hi/ A. Let 
us nevertheless proceed and calculate the corresponding 
dynamic structure factor, which can be obtained by re- 
placing the term hi/ A = hi/(xQ + y 2 ) on the right-hand 
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9 = 1 
Po = 0.2 
p = 0.04 




X — Xq 



FIG. 9: (Color online) Graph of the dynamic structure factor 
S(lo, q) as a function of x — xn = (u) — vq) /(uf?) for fixed q = 
0.08/cf. The line-shape has been calculated from Eqs. (|5.18l 
I5.20|l and ()5.22[) . The distance between the local maxima is 
proportional to w 9 tx q 2 /m. 



of the form f q = f + ±/^g 2 + 0(q 4 ) with fg ^ we do 
not encounter any mass-shell singularities. In fact, we 
believe that even for sharp momentum-transfer cutoff, 
fq = /o©(<7o — ?)i our perturbative result (|4.ip does not 
suffer from mass-shell singularities as long as we do not 
rely on the approximation A discussed in Sec. IIVB1 in 
other words, the mass-shell singularity hi / A in Eq. (|5. 18|) 
is an artefact of the sharp momentum-transfer cutoff in 
combination with our neglect of curvature corrections to 
the free polarization in loop integrations. While we are 
not able to evaluate Eqs. (|4. 14114.17)) analytically with- 
out relying on approximation A, we shall in this sec- 
tion abandon the sharp momentum-transfer cutoff and 
assume that the interaction f q can be expanded for small 
q as in Eq. (|1.15p . Later we shall argue that as long as we 
rely on approximation A, our result for S(u>, q) can only 
be trusted for q > q c = l/(m|/o'|), see Eq. ([Tig]) . But 
if /q is sufficiently large, then there exists a parametri- 
cally large regime q c -c q <C kp of wave- vectors where 
our calculation is valid. 



side of Eq. IpTTgjt by 
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(5.22) 



The finite imaginary part w g in this expression is a rough 
estimate of the modification of the spectral line-shape 
due to the terms which have been neglected by making 
the approximation A discussed in Sec. IIV Bl The typ- 
ical form of the dynamic structure factor in the regime 
p < po implied by Eqs. (|5. 18(15. 20|) and (|5.22p is shown in 
Fig. [9l Obviously, within our approximation the dynamic 
structure factor does not exhibit any threshold singular- 
ities, which according to Refs. [1Q are a generic feature 
of the dynamic structure factor of Luttinger liquids. It 
turns out that the absence of threshold singularities in 
Fig. [9] is an artefact of the rather simple rcgularization 
prescription (|5.22p of the unphysical mass-shell singu- 
larity in Eq. (|5. 18[) . In the following section we shall 
show how to recover the threshold singularities within 
our functional bosonization approach. 



A. Imaginary part of II, 1 (uj,q) 



Let us first calculate the imaginary part of the di- 
mcnsionless inverse polarization II~ 1 (a; + i0,p) given 
in Eq. (|4.44p assuming for simplicity p > 0. From 
Eqs. P~4"5)) and P~4"?j) we obtain 



Imli,, 1 (x + iO, p) = Im/(x + iO, p) 
1 



dp p Im 



J(x + i0,p,p') + J(—x - i0,p,p') 



(6.1) 



VI. INTERACTION WITH REGULAR 
MOMENTUM DEPENDENCE 

In this section we shall show that for a more realistic 
interaction whose Fourier transform is for small momenta 



In order to calculate the imaginary part of J(x+i0,p,p'), 
we first perform a partial fraction decomposition of 
Eq. (|4.47p . then carry out the analytic continuation to 
the real frequency axis iy — > x + iO, and finally take the 
imaginary part using Im[a — x — i0] _1 = irS(a — x). After 
some lengthy algebra we obtain 
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lmJ(x + iO,p,p') = — 
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(6.2) 



where we have defined i p + p ' = sign(p+p')x p + p / . In order 
to perform the p'-integration in Eq. (|6.ip , we use the fact 
that by assumption both p and p' are small compared 
with unity so that we may expand x p to first order in p 2 , 



x p = v 71 + 9 P 



x 



where from Eq. (|4.11[) . 
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Note that for small p c the coefficient x'q is large compared 
with unity. The (^-functions in Eq. (|6.2[) can then be 
approximated by 
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(6.6) 



In Eq. (|6.5[) we have expanded the argument of the 5- 
function to linear order in p and p' , assuming that both 
dimensionless momenta are small. On the other hand, 
due to the cancellation of the leading term in the dif- 
ference Xp' — x p i + p in the (5-function of Eq. (|6.6[) . the 
corresponding expansion has to be carried out to cubic 
order in the momenta. The integration in Eq. (|6.2p can 
now be carried out analytically and we obtain for small 
P >0, 
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FIG. 10: (Color online) Graph of the functions Ci(u) and 
C R {u) defined in Eqs. fT9)) and (|6TT3)| . The dotted lines in- 
dicate asymptotic limits. 



where 



7p 



3p 2 



8irx Pc 

and the function Ci (u) is given by 
d(u) = e(u)Q(l - u) 



(6.8) 



(6.9) 



Note that the coefficient hi = (1 + 3xq) 2 / (12xq) on the 
right-hand side of Eq. (|6.7j) has also appeared for sharp 
momentum-transfer cutoff [see Eq. (|5.19c[) ] in form of 
the residue of the mass-shell singularity h\j A in our 
expression (|5 . 1 8(1 for the irreducible polarization. A 
graph of Cj(u) is shown as the dashed line in Fig. [TO] 
Mathematically, the square-root singularity of Ci(u) for 



1 originates from the special point x- 



7 P where 



the argument of the Dirac (5-function on the right-hand 
side of Eq. ()6.5j) has a double root. We believe that the 
divergence of Ci(u) for u — > 1 is unphysical and indicates 
that the approximations leading to Eq. (|6.5[) are not suf- 
ficient in this regime. Hence, within our approximations 
we can only obtain reliable results for the spectral line- 
shape as long as the ratio (x — x P )/7 P is not too close to 
unity. 
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B. Real part of IR 1 (uJ,q) 

For p c -C 1 and p <C 1 we can obtain the contribu- 
tion from RcI(x + iO,p) analytically from Eqs. (|4.45p and 
(|4.47p using the fact that among the corrections of order 
p 2 only terms proportional to p 2 jp c should be retained. 
We obtain for x > and p > 0, 



Rel(x + i0,p) = h-x 2 I 2 



7rp c /iisigii/o'C,R 



7 P 



(6.10) 



where 
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(6.11) 
(6.12) 
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and the function Cr(u) is given by 
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6(1 -u) In 
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A graph of Cr(u) is shown in Fig. [TO] (solid line). 
Note that Cr(u) and Cj{u) can be written as Cr(u) = 
ReC(u-HO) and Ci{u) = ImC(u-MO), where the complex 
function C{z) is 



C(z) 
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The real part of our dimcnsionlcss inverse polarization 
can be written as 

Ren^x-HO,^) = Z x -Z 2 x 2 

+ npchtm^fdcJ^Z^), (6.15) 



with 



Zi = l + h 



'H, 



Z 7 



(6.16) 



By assumption, the bare interaction f q is negligibly 
small for momentum-transfers exceeding qo <C kp 7 so 
that the integrals I\, I 2 and Ih are proportional to 
Pa = [<7o/( 2fc F)] 2 < 1 and hence Zi = 1 + 0(pl). Keep- 
ing in mind the self-consistent definition ()4.8|) of xo, we 
finally obtain for positive x and p, 

g p + Refl- 1 (x + iO,p) = Z 2 [x 2 p -x 2 + R(x,p)], (6.17) 

where 

= ^signtfcJ^Y (6.18) 
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FIG. 11: (Color online) Graph of the dynamic structure fac- 
tor S(u>,q) given in Eq. (|6.21[) as a function of x — x v for 
p = 0.04 = 25p c and g = 1. For simplicity we have set 
Zj » 1, which is accurate for po <C 1. For p 3> p c most 
of the spectral weight is carried by the main shoulder whose 
lower edge bounded by a threshold singularity. The 

width of the main shoulder on the x axis scales as 7 P oc p 2 /p c - 
Recall that x — u}/(vpq), so that the corresponding width on 
the frequency axis scales as j q — VFqjp oc q 3 /(mq c ). For 
p p c the small "satellite peak" emerging above the upper 
edge of the main shoulder carries negligible spectral weight 
and is probably an artefact of our approximations. 



C. Spectral line-shape of S(uJ,q) 

To discuss the line-shape of the dynamic structure fac- 
tor, it is convenient to introduce also the imaginary part 
of the effective self-energy via 



Imfi- 1 ^ + z0,p) = -Z 2 T(x,p), 



(6.19) 



or explicitly 

T(x,p) 
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The dynamic structure factor can then we written as 



S(u,q) 
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The resulting line-shape for p p c is shown in Fig. 1111 
Obviously, S(u>, q) exhibits a threshold singularity at x = 
x p , corresponding to the threshold frequency 



v F qx p = vq + 



signJo q 3 
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Moreover, most of the spectral weight is smeared out over 
the interval < x — x p < 7 P , or equivalently lo^ < lu < 



7 g , where the energy scale 7 g is defined by 



lq = v F qj p 



8ttxq 2mq c 



(6.23) 
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FIG. 12: (Color online) Solid line: dimensionless ZS damp- 
ing 7 P = 7(j/(«F?) denned in Eq. (|6.23|l as a function of 
p/p*. Dashed line: estimate of the width w p = w q /(vFq) = 
(Z w /V3)p of the ZS resonance given in Eq. (|5 . 20 fl . 



The energy "f q can be identified with the width of the 
ZS resonance on the frequency axis. The crucial point 
is now that for q » q c Eq. (|6.23p is much larger than 
the estimated broadening w q cx q 2 /m of the ZS reso- 
nance due to the terms which we have neglected by mak- 
ing the approximation A discussed in Sec. IIVBI (which 
amounts to ignoring non-linear terms in the energy dis- 
persion in bosonic loop integrations). Our approxima- 
tion A is therefore only justified in the regime where the 
broadening j q due to the q-dependence of the interac- 
tion f q is large compared with the broadening w q due 
to the non-linear energy dispersion in bosonic loop inte- 
grations. We thus conclude that the calculations in this 
section are only valid as long as j q > w q . A comparison 
of 7 g and w 9 is shown in Fig. 1121 Obviously, the condi- 
tion w q — 7 9 defines a characteristic crossover scale q* 
where the g-dependence of the width of the ZS resonance 
changes from q 2 to q 3 . Using Eqs. (|5.20j) and (|6.23p we 
obtain the following estimate for the crossover momen- 
tum scale, 



8ttZ w xq 
3V3 
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which has the same order of magnitude as q c = 
l/(m|/o'|). We conclude that the results for S(u>,q) pre- 
sented in this section are only valid for q > q* , and hence 
do not describe the asymptotic q — > regime. But the 
scale q r can be quite small for some interactions. For 
example, if the interaction f q can be approximated by a 
Lorentzian (|1.13[) with screening wave- vector qo <C hp, 
then q c = gg/(2m/o) is quadratic in qo. For long-range 
interactions the regime q* < q -C qo where our calcu- 
lation is valid can therefore be quite large and physi- 
cally more relevant than the asymptotic long-wavelength 
regime g « g,. 

The small "satellite peak" slightly above the main 
shoulder in Fig. [11] is probably an artefact of our approx- 



imations, in particular of approximation A discussed in 
Sec. IIV Bl It is easy to show that the satellite peak is lo- 
cated a distance Sx cx p\/p 2 <C 7 P above the upper edge 
Xp + jp of the main shoulder and its width is proportional 
to p 2 7 P cx p^/pc <C Sx <C 7 p . Note that in the regime 
q 3> q c where our calculation is valid the threshold sin- 
gularity is located at tu~ ~ vq — 7 9 (up to corrections of 
the order q 2 /m <C 7 g ), while the energy scale of the satel- 
lite peak is vq + 0(q 2 /m). However, as discussed after 
Eq. ()6.9|) . in the regime \ (x — x p ) / % — 1| <C 1 our approx- 
imation A is not reliable, so that the detailed line-shape 
in the vicinity of the satellite peak is probably incorrect. 
Fortunately, for p 3> p c the satellite peak carries negligi- 
ble weight, so that our calculation reproduces the main 
features of the spectral line-shape. We speculate that a 
more accurate evaluation of our self-consistency equation 
for n*(<j, q) derived in Sec. IIV Ai which does not rely on 
approximation A in Sec. HVBi will generate additional 
weight in the dip between the upper edge of the main 
shoulder and the satellite peak, resulting in a single local 
maximum at the upper edge of the main shoulder. The 
spectral line-shape looks then qualitatively similar to the 
line-shape proposed in Refs. 

Let us next consider the line-shape in the vicinity of the 
threshold singularity x — > x p . For < {x — x p )/"/ p <C 1 
we may approximate 
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In the last line we have approximated Z 2 ~ 1. From the 
above discussion it is clear that this expression can only 
be trusted for p > p* . A graph of function of p/p* 

is shown in Fig.[T3J In the regime \rj v \ ln[47 P /(x — x p )] 
1, which is equivalent with 



< x - x p < 47 P exp [— l/|r/ p |] 



(6.28) 



the dynamic structure factor can thus be approximated 
by 
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FIG. 13: (Color online) Solid line: graph ol r\ v denned in 
Eq. (|6.27p as a lunction ol p/p* for /,'/ < 0. The dashed line 
is the weak coupling result r\ v w 1/2 — p/(47rp c ) obtained by 
Pustilnik et al. in Rel. [J. The dashed dotted curve is a 
simple parabolic interpolation. 



According to Pustilnik et alA, the logarithmic singularity 
can be re-summed to all orders, so that it is transformed 
into an algebraic one. Assuming that this is indeed cor- 
rect, we can replace 



R(x,p) 
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For x 
as 



— > 2xo(x — .T p ) 
x p the dynamic structure factor then diverges 
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2.t Z 2 (47p) 2 " P [a - x p ]^ 
with the threshold exponent 



(6.31) 



calculated exactly^ the momentum-dependence of ?y p 
looks different from our result for the FSM. For exam- 
ple, in the Calogero-Sutherland model r] p is independent 
of p, see Ref. However, the Fourier transform f q of 
the interaction in the Calogero-Sutherland model van- 
ishes for q = 0, while in the intcgrablc XXZ-chain con- 
sidered in Refs. f^l MIl^lTlT ] the effective interaction of the 
equivalent one-dimensional fermion system involves also 
momentum-transfers of the order of Moreover, in the 
XXZ-chain there exists no crossover scale q c satisfying 
q c = (m|/o |) _1 <C kp, so that the intermediate regime 
q c <C q -C kp where j q oc q 3 /q c simply does not exist. 
The existence of such an intermediate regime seems to be 
a special feature of the FSM considered here, where f q 
involves only small momentum-transfers and has a finite 
limit for q = 0. 

Within our perturbative approach we cannot justify 
the re-summation procedure (|6.30|) . Possibly a care- 
ful analysis of the functional rcnormalization group flow 
equation for the irreducible polarization discussed in ap- 
pcndix|D]will shed some light onto this difficult problem. 
This seems to require extensive numerics, which is be- 
yond the scope of this work. Note that for /g > the 
exponent 7] p in Eq. (|6.27p is negative, so that the sin- 
gularity in Eq. (|6.31[) is not integrable and exact sum 
rulesi£ cannot be satisfied. In contrast, the original log- 
arithmic singularity in Eq. (|6.29|) is integrable (the inte- 
gral J Q dt/[thi 2 t] is finite), so that at least for /q > 
the logarithm found in perturbation theory cannot be 
exponentiated. On the other hand, an interaction with 
/o > seems to be unphysical and does not describe a 
stable Luttinger liquid^. 

Finally, consider the tails of the spectral function. For 
x 3> x p we obtain from Eqs. (|6.20|) and (|6.21l) . 
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Note that for /q < and p <C 1 the weak coupling esti- 
mate for fx p given by Pustilnik et alA is in our notation 
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Inserting our result (|6.8[) for 7 P we obtain 
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As shown in Fig. [T3l this is consistent with a smooth 
crossover to our result (|6.27p at p/p* = 0(1). Qualita- 
tively, we expect that the behavior of rj p in the crossover 
regime resembles the dashed-dotted interpolation curve 
in Fig. [T31 Note that T) v < 1/2 for all p, so that p, p > 0. 
For some intcgrablc models where rj p has recently been 



S(u,q) 
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(6.37) 



in agreement with Refs. Note that the tail of 

S(uj,q) is determined by the first term on the right-hand 
side of the damping function T(x,p) given in Eq. (|6.20|) . 
whereas the regime close to the ZS resonance is deter- 
mined by the second term involving the complex func- 
tion C(z). This is the reason why the spectral line-shape 
close to the ZS resonance cannot be simply obtained via 
extrapolation from the tails assuming a Lorentzian line- 
shape. 
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VII. SUMMARY AND CONCLUSIONS 

In this work we have used a functional bosonization ap- 
proach to calculate the dynamic structure factor S(u,q) 
of a generalized Tomonaga model (which we have called 
forward scattering model) , consisting of spinless fermions 
in one dimension with quadratic energy dispersion and 
an effective density-density interaction involving only 
momentum-transfers which are small compared with kp. 
We have derived in Sec. lIVI a self-consistency equation for 
the irreducible polarization II* (u>, q) which does not suf- 
fer from the mass-shell singularities encountered in other 
perturbative approaches. Although for the explicit eval- 
uation of S(to,q) we had to make some drastic approx- 
imations (in particular, in bosonic loop integrations we 
have neglected curvature corrections to the free polar- 
ization, see approximation A discussed in Sec. IIV Bp we 
have found a regime of wave- vectors q c <C q <C kp where 
an explicit analytic calculation of the spectral line-shape 
is possible. The crossover scale q c = l/(m|/o'|) is deter- 
mined by the second derivative /q' of the Fourier trans- 
form of interaction at q = 0. For interactions whose 
Fourier transform can be approximated by a Lorentzian 
with screening wave- vector qa <C kp , the crossover scale 
q c is proportional to q 2 , so that the regime q c <C q <C kp 
is quite large and can be experimentally more relevant 
than the asymptotic long- wavelength regime q <C q c - We 
have shown that for q c <C q <C kp the width of the ZS 
resonance on the frequency axis scales as 7 g cx q 3 / (mq c ). 
Our result is consistent with a smooth crossover at q f=a q c 
to the asymptotic long- wavelength result 7 g oc q 2 jvn ob- 
tained by other authors^^. The spectral line-shape is 
non-Lorentzian, with a main hump whose low-energy side 
is bounded by a threshold singularity at w = uj q = 
vq — 7 9 , a small local maximum around to vq, and a 
high-frequency tail which scales as q A /oj 2 . For w — > uj^+Q 
the threshold singularity is within our approximation log- 
arithmic, S(lo, q) cx [(lu — u> q ) ln 2 (cj — lo^)]^ 1 . Assuming 
that higher orders in perturbation theory exponentiate 
the logarithm, we obtain an algebraic threshold singu- 
larity with exponent [x q = 1 — 2rj q and rj q oc q 2 /q 2 for 
q > q c - 

Finally, let us point out a number of open problems: 

1. It is by now established that, at least in inte- 
grable models, S(ui, q) indeed exhibits algebraic thresh- 
old singularities^ ^ 10 ' 11 . However, for generic non- 
integrable models there is no proof that the logarithmic 
singularities generated in higher orders of perturbation 
theory indeed conspire to transform the logarithm en- 
countered at the first order into an algebraic singularity, 
as suggested in Ref. Q. This would require a thorough 
analysis of the higher order terms in perturbative expan- 
sion, which so far has not been performed. 

2. For the explicit evaluation of the self-consistency 
equation for the irreducible polarization II* {uj, q) derived 
in Sec. IIV Al we had to rely in this work on approxima- 



tion A discussed in Sec. IIV Bl We have argued that this 
approximation is not sufficient to calculate the dynamic 
structure factor for q < q c , because it neglects the dom- 
inant damping mechanism in this regime. Moreover, for 
sharp momentum-transfer cutoff our approximation A 
breaks down for frequencies in the vicinity of the mass- 
shell singularity. It would be interesting to evaluate the 
self-consistency equation for the irreducible polarization 
Tl*(uj,q) derived in Sec. IIV Al without relying on approx- 
imation A. We believe that in this case our functional 
bosonization result for S(uj, q) does not exhibit any mass- 
shell singularities even for sharp cutoff. The explicit eval- 
uation of the relevant integrals is quite challenging and 
probably requires considerable numerical effort (includ- 
ing a numerical analytic continuation), which is beyond 
the scope of this work. 

3. By assumption, the interaction of the FSM con- 
sidered in this work is dominated by small momentum- 
transfers q <^ kp. On the other hand, the Fourier trans- 
form of the effective interaction in the Jordan- Wigncr 
transformed XXZ-chain studied in Refs. UGlIllI has also 
components involving momentum-transfers of the order 
of kp. It should be interesting to investigate more thor- 
oughly how the dynamic structure factor depends on the 
properties of the interaction. Unfortunately, the FSM 
discussed in this work is not integrable and there seems to 
be no integrable model with quadratic energy dispersion 
where the interaction involves only small momentum- 
transfers and has a finite limit for q — > 0. 

4. In appendix [D] we present a functional renormal- 
ization group equation [see Eq. (|D14[) ] for the irreducible 
polarization which goes beyond the self-consistent per- 
turbation theory based on functional bosonization used 
here. A thorough analysis of Eq. (|E)14|) using numerical 
methods still remains to be done. Possibly, this equation 
will be a good starting point for addressing some of the 
open problems mentioned above. 
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APPENDIX A: FERMION LOOPS FOR 
QUADRATIC DISPERSION IN ONE DIMENSION 

In the functional bosonization approach the vertices 
of the interaction part S- m t[5(j)] of the bosonized action 
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((3~2"0|) are 

I^Ql, . . . , Q n ) = l n {n - 1)! 4 n) (-«l. • • ■ . ~Qn), 

(Al) 

where the symmetrized closed fcrmion loops are defined 
by 

L s (Qx, . . . , Q n ) 



1 E / G (K)G (K-Q P{1) 

% ' P(l,...,n) K 



n-1 



xGo(A' - Q P(1) - Q P(2) ) ■ • • G (K - QpU)) ■ ( A2 ) 

3=1 

Here the sum is over all permutations P(l, ...,n) of 
1, . . . , n, and the fermionic Green functions Go (if) should 
be calculated within the self-consistent Hartree approx- 
imation, see Eq. (|3.16p . For fermions with quadratic 
energy dispersion in one dimension, the symmetrized 
fermion loops (|A2[) can be calculated exactly. Neumayr 
and Metzne r 28 i 29 (see also Ref. [3(|) have derived reduc- 
tion formulas for quadratic dispersion in D dimensions 
which allow to express the non-symmetrized loops 

« n 

L^\Q 1 ,...,Q n ) = J^llGoiK-Qt) 



G (K-Q 1 )Go(K-Q 2 )---G (K 



0, (A3) 



K 



for n > D + 1 in terms of linear combinations of 
the more elementary loop L^ D+1 \Qi, . . . , Qd+i)- In 
particular, in D = 1 the non-symmetrized loops 
L^(Qi, . . . , Q n ) with n > 2 can be expressed in terms 
of the two-loop Z( 2 )(0, -Q) = L ( g \~Q,Q) = -n (Q). 
Given explicit expressions for the non-symmetrized loops 
L^ n '(Qi, . . . , Q n ) we may construct the corresponding 
symmetrized loops L^(Qi, . . . ,Q n ) by shifting the la- 
bels, 



Qi 


= o, 


Q~2 


= Qi, 


Q:< 


= Qi + C 


Qn 


n-1 

= 

3 = 1 



(A4) 

so that Qi+i — Qi = Qi, and defining 

L ( - n \Q 1 ,...,Q n ) = L ( - n \Q 1 ,...,Q n ). (A5) 
Then the symmetrized loops are 

( " w n "-^ E ^Wpw,-.^-))- 

(A6) 



£g (Ql, ■ • ■ , Qn) 



P(l,...,n) 



In one dimension, the reduction formula for the non- 
symmetrized loop L^(Qi, . . . , Q n ) given by Neumayr 
and Metzner— can be obtained using a straight-forward 
partial fraction decomposition. Performing the frequency 
integration in Eq. (|A3[) and introducing the notation 
Qi = (iuii, qi) we obtain 



£ (n) (Qi 



o = E 



k [7 



dk "I r 1 

27T 

3#« 



(A7) 



where 



£lij(k) = i(bJi - (Jj) + & - Ck+q.-q 



'i j ■ 



(A8) 



and & = ^ + /oPo - M = ( fc2 - &|)/(2m). Defining 



A'; 



9j Qi , ■ 
'■ vm 



F i 
UJi — UJi 



2 % - 

we may alternatively write Eq. (|A7|) as 



l(«)(q 1 ,...,q„)=^/ 



(A9) 



~ ./-fc F 27r ^ {Qj ~ Qi)(k - hj) ' 

(A10) 

We can now perform another partial fraction expansion 
to obtain 



L (n) (Qi,...,Q„) = J2 



• ll 



qj - * 



:f dk 1 



(All) 



with 



„ fa - ii)(<ii - Qj) , y- — \ . ■ / — -,qi-q 



In the special case n = 2 this yields 



L (2) (Oi,Qa) 



7r(gi - g 2 ) 



In 



A;f — fci 



% - ft 
(A12) 



(A13) 



In order to give an explicit formula for the function 
L^ n \Qi, . . • , Qn) defined in Eq. (|A5[> . which depends on 
the external momenta and frequencies Qi = (iu)i,qi), we 
introduce the notation 



a- - a a - 1 ^'=J 9 ' ' * > J 

I - E/ =i qu j > i 



(A14) 



and similarly for un 



ujj. These quantities fulfill 



qij = qn + qij and (Vij = ll>u + u>ij , such that H^i can be 
reexpressed as 



H 



1 r 

q>j 



i(vaqij - quuiij) 



quqijqij 



2m 



(A15) 
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which is manifestly symmetric under exchange of i and 
j, i.e., Hju = H ijt . This yields 



= - 



i<j 



II 



i = l 



n (Qi 



(A16) 

with = (iujjj , qjj) . This result is equivalent with 
Eq. (19) of Ref. [20| . Finally, in order to obtain the sym- 
metrized loops Lg (Qx, . . . ,Q n ) in Eq. (|A6|) . an addi- 
tional summation over the n\ permutations is necessary. 
Evidently, the resulting expressions are rather compli- 
cated. In the following two appendices we shall therefore 
discuss the symmetrized three-loop and the symmetrized 
four-loop separately. However, without explicitly evalu- 
ating the loops the following two general properties can 
be established: 

1. The symmetrized n-loops Lg (Qi, . . . , Q n ) are fi- 
nite for all values of their arguments 2 ^. This guar- 
antees that in the perturbative expansion of the 
irreducible polarization Tl*(Q) in powers of the 
RPA interaction no infrared singularities are en- 
countered. 

2. In the limit 1/m — > the symmetrized n- 
loop is proportional to (l/m)"~ 2 . More pre- 
cisely, the dimcnsionlcss symmetrized n-loops 



Lq (Qi, ■ ■ ■ j Qn), defined via 
(n-l)!4" ) (Q 1 ,...,Q„ 



(A17) 



have finite limits for 1/m — * 0. For large m the 
vertices r^^i, . . . , Q n ) in the interaction part 
Sint[<^] of our effective action (|3.20p are therefore 
proportional to increasing powers of the small pa- 
rameter 1/m, which justifies the perturbative treat- 
ment of these vertices. 



APPENDIX B: SYMMETRIZED THREE-LOOP 

The explicit expression for the symmetrized three-loop 
can be written as 



(3) 
J s 

= -Re 



Ly (iwi, qi;iu 2 , qi\ -iui - iu 2 ,-qi - q 2 ) 
1 



iLUiq 2 - iuj 2 qi - qiqi^jr 



qiUo(iui,qi) + q 2 Tl (iLU 2l q 2 ) 
~(qi + q 2 )Tl Q (iuJi + ioj 2 , qi + q 2 ) 



(Bl) 



Introducing again the variables iy± = iu>x/(vFqi), Pi = 
qi/(2kp) (and similarly for iy 2 and p 2 ) and the dimen- 
sionlcss function ILo(iy,p) = Uo(iu>, q) [sec Eqs. (|2.11[) 



and (|2.12[) ] , we may write the symmetrized three- loop in 
the dimensionless form (|A17|) . 



2L i §\iui,qi 1 iuj 2 ,q 2 - 1 —iwi - iu 2 ,-qi - q 2 ) 

= — 2 L s '(iyi,pi]iy 2 ,p 2 ) , 
mvp 



(B2) 



with 



Ls{iyi,pi;iy 2 ,p 2 ) 



(yi - V2) 2 + {pi +p 2 ) 2 



l - . , l - . 

— n (zj/i,pi) H Il a (iy 2 ,p 2 ) 

s 2 Si 



1 



si s 2 



fio(%si + iy2s 2l pi +p 2 ) 



, (B3) 



where we have defined 
Pi r 



si 



with 



P1+P2 r + 1 



s 2 



P2 



Pl+P 2 T + l 



r = — 
Pi' 

For later convenience we also define 
Pi r p 2 



, (B4) 



(B5) 



n = 



P1-P2 r - 1 



P2-P1 r — 1 



(B6) 



Note that by construction si + s 2 = T\ + r 2 = 1. 

At the first sight it seems that the symmetrized three- 
loop diverges for \pi/p 2 \ — » or \p 2 /pi\ — > 0. More- 
over, the prefactor in Eq. (|B1[) diverges in the special 
limit pi — > p 2 and f/j — > y 2 . It turns out, however, 
that all divergencies cancel and the symmetrized three- 
loop is everywhere of the order of unity. This non-trivial 
cancellation cannot be obtained by power-counting and 
can be viewed to be a consequence of the asymptotic 
Ward-identity associated with the separate conservation 
of left-and right-moving particles for linearized energy 
dispcrsio n 32 i 33 . We shall show shortly that a similar can- 
cellation protects also the symmetrized four-loop from 
divergencies. The symmetrization of the loops is crucial 
to cancel the divergencies. In diagrammatic language, 
the symmetrization properly takes vertex and self-energy 
corrections into account. 

The limiting behavior of the function 



Ly(iyi,pi;iy 2 ,p 2 ) for p x 



and p 2 



is not 
Using 



unique but depends on the ratio r — pi/p 2 . 
Eq. (|2.16[) we obtain after some algebra, 



lim L i p(iy 1 ,p 1 ;iy 2 ,p 2 ) = L s %{iy 1 ,iy 2 ,r) 
-»0,pi /P2— r 



with 



L%{iyi,iy 2 ,r) = 

1 - 3/12/2 - {m + j/gXfij/i + s 2 y 2 ) 
[I + y 2 ][l + yl][l + {siyi + s 2 y 2 y] 



(B7) 



(B8) 



28 



1)2 



AO 




in 



FIG. 14: (Color online) Graph of the function L s o {iyi, iy2,r) 
given in Eq. (|B8|I for ya = 1 as a function of j/i and r = Pi/p2- 



which is manifestly finite for all values of its arguments. 
A graph of the function L s (iyi, iy2, r) is shown in 
Fig. [14] Note that a finite limit of the dimensionlcss 

function L s ' {iyi,pi;iy2,P2) for small momenta does not 
contradict the loop cancellation theore m 25 i 27 i 28 ' 29 i 32 ' 33 , 
because according to Eq. (|B2[) the physical symmetrized 

three-loop Ly(iui, qr, iu 2 , 92; -i^i - i^2,—Qi - 92) in- 
volves an extra factor of 1/m, so that it vanishes for 
1/m -> 0. 



APPENDIX C: SYMMETRIZED FOUR-LOOP 

The symmetrized four-loop is more complicated than 
the three-loop. However, the four-loop determines the 
correction to the irreducible polarization to first order in 
the RPA interaction, so that we need it for our calcula- 
tion. Actually, we need the four-loop only for the special 
arguments Q3 = —Q\ and Q4 = —Q%. It is useful to 
introduce the notation 



y± 
p± 

and the complex functions 

C±{iy-,pi,pz) = 



yi ± 2/2 , 

Pl ± P2 , 



P1P2W- 



P±\ 



(Cla) 
(Clb) 



(C2) 



W(iy,p) 
We also need 



1 

2^ 



1 



1 



iy + 1 + p iy + 1 — p 



ReW(iy,p) = 
lmW(iy,p) = 



y 2 -i+P 2 

[y 2 + (l+p) 2 ][y 2 + (l-p) 2 ] 

2y 

]y* + (i+p) 2 W + {i-p) 2 ] 



(C3) 

(C4) 
(C5) 



The functions C±(iy-,pi,p2) are singular for pi — > 0, 
while W(iy,p) has a finite limit for p — > 0, 



]im W(iy,p) = - 

p~*0 (1 + iy) 2 - 



(C6) 



Using our general result (|A16[) for the non-symmetrized 
n-loops L^ n \Qi, . . . , Q n ) and performing the sum (|A6|) 
over all permutations of the external labels, we obtain 
for the dimensionlcss symmetrized four-loop (as defined 
in Eq. (|A17|) for n = 4) for the special combination of 
external labels needed in Eq. (|4.14[) . 



6Ly(iuji,qi; —iui, -qi;iuj 2 , 92; -i^2, -92) 

„ v o ?(4), 



(mvjr) 



2 \2 L S (%)2>i;*2fe,i>2) J 



(C7) 



with 



Ls Hm,piim,P2) = 

+^-r c [p+c 2 + P _c 2 + 2 Pl c;c_] n (%,pi) 

+^Re [p + C 2 -p_C 2 - 2p 2 C + C_] n (zy 2 ,P2) 

-p 2 _[ReC+] 2 n (i?/iSi +iy 2 s 2 ,P+) 
-p 2 _[ReC-] 2 IL (iyin + 12/2^2, P-) 

+ ilm[C+ - C-JIm^Ciy!^!) - W(iy 2 ,p 2 )] 
-Re[W(iy 1 ,p 1 )W(iy 2 ,P2)] , (C8) 



where we have written C± = C±(iy-,pi,p 2 ). After some 
algebra Eq. (|C8|) can be cast into the form 
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Ls'(iyi,pr,iy2,P2) = 2 2T5 - 1 — 2m 2 I — 2I2 n (iyi,pi) 

P2W- +P+} 2 [y- + P A]2 



M + ^y- + ^{pl-p\?y- + (pl-pj) 3 ^ ,. s 

P 2 [y 2 +P 2 g 2 [ 2 / 2 +P 2 J 2 n °(M) 

n (iyisi + iy2S2,p+) tlo(iyin + iy 2 r 2 ,pJ) 



s 2 s 2 ,[y 2 +P+] 2 rl4[y 2 _ + p 2 _} 2 

2y_lm[W(iyx,px) - W(iy 2 ,p 2 )] 



[y-+p\][y-+p-] 
1 



Re[W{iy 1 ,p 1 )W(iy 2 ,P2)] ■ (C9) 



Naive power counting would suggest that this expres- 
sion is singular for 2/1 — ► 2/2 or \p\\ — ► |jp2 1 , or if p\jp 2 
approaches either zero or infinity. However, similar to 
the symmetrized three-loop, all singularities cancel in 
Eq. (|C9[) . so that the symmetrized four-loop remains fi- 
nite and of the order of unity for all values of its argu- 
ments. 



For simplicity, consider again the limit p\ — > and 
P2 — > with constant r = p\/p2- Then 

lim L { g\iy 1 ,p 1 ,iy 2 ,P2) = Ly {iy 1 ,iy 2 ,r) , 

Pi— ►0,pi/P2=r 

(CIO) 

with 



L { sl(iyi,iy2,r) 



1 



(yi - 2/2)" 



2(r 2 +3) 2(r- 2 +3) 



i + yl 



i + yl 



s\sl[l + (siyi + S2J/2) 2 ] rfr^l + (nyi + r 2 y 2 ) 2 



4y 2 [1 - 2 yi y 2 - 2/12/2(2/1 + y\ + ym)} 



[i + y 2 ] 2 [i + 2/ 2 2 ] 2 



4yiy2-(l-y 2 )(l-2;2 2 ) 



[i + y 2 ] 2 [i + y 2 2 ] 



212 



(Cll) 



The important point is now that the singular prefactor l/(j/i — y 2 ) 4 m Eq. (|C11[) is compensated by a factor of 
(2/1 — y 2 ) 4 arising from the sum of the five terms in the square braces. In fact, we can explicitly cancel this singularity 
by combining these terms differently, 



L%{iyi,iy2,r) 



[i-y 2 ][i-y 2 2 ] 



[1 + y 2 ] 2 [l + y 2 } 2 [1 + y 2 ] 2 [l + y 2 2 ] 2 [l + (s m + s 2 y 2 ) 2 ][l + (nyi + r 2 y 2 ) 2 ] \ 
2r 2 s 2 [l + y 2 ][-l + 3y 2 + y 2 + 82/12/2 - 3yjy 2 } + 2r 2 s 2 [l + y 2 ][-l + 3y 2 + y 2 + 8 yi y 2 - 3y 2 y 2 } 



Hrlsi, + tiai) [1 + yf][l + l^][-6 + y{ + 2/2 + 4yiy 2 ] + 8[1 - 2y lV2 - Syfyi] 



(C12) 



Note that the coefficients fx, r 2 , si and S2 are not independent but can be expressed in terms of a single parameter 
r = pi/p2, as given in Eqs. (|B4[) and (|B6|) . Introducing the notation 



h = siri 



Pi 



p\-p\ r ' 2 - 1 
so that ti + £2 = 1 ; wc may alternatively write 

[l-2/ 2 ][l-2/ 2 2 ] 



, *2 = s 2 r 2 



pI 



9 9 1 5 

P 2 - PI - 1 



(C13) 



L { s\{iyi,iy2,r) 



[1 + y 2 ] 2 [l + 2/f] 2 [1 + ?y 2 ] 2 [l + y|] 2 [l + (s iyi + s 2 y 2 ) 2 }[l + (ru/i + r 2 y 2 ) 2 } \ 
-1 + 62/12/2 +£1*2(2/1 - 2/2) 2 [j/ 2 + y\ + 62/12/2] 

+2(ti2/i + t22/2) 2 2/i2/2(4 - 2/12/2) + 2(ii2/i + t 2 y 2 ) - t 2 y 2 )(yl - y 2 ) + + *2j/i) 



(C14) 
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1/2 = 1 




r 



FIG. 15: (Color online) Graph of the function Lg (iyi, iy2, r) 
given in Eq. (|C11|I for y2 = 1 as a function of y\ and r — pi jpi ■ 



A graph of the function L s (iyi, iy2: r) is shown in 
Fig.HU 



APPENDIX D: FUNCTIONAL 
RENORMALIZATION GROUP EQUATION FOR 
THE IRREDUCIBLE POLARIZATION 

In this appendix we shall derive a formally exact func- 
tional rcnormalization group (FRG) equation for the ir- 
reducible polarization which can be used to generate a 
perturbative expansion of the polarization in powers of 
bosonic loops with either the RPA interaction or the 
true effective interaction as propagators. We use the 
momentum-transfer cutoff scheme proposed in Ref. [36j , 
where only the free bosonic part Sq[4>] is regularized by 
suppressing bosonic fluctuations with momenta q smaller 



than a certain cutoff A. One possibility is to introduce 
the cutoff as a multiplicative 0-functionsS by replacing 
in Eq. ^ 

f q G(q > \q\ > A)/„ (Dl) 

where 0(A) = 1 if the logical expression X is true, and 
0(A) = if X is wrong. Alternatively, we may insert an 
additive cutoff R\(q) in the inverse propagator—, 

f-'^f-'+RAiq), (D2) 

where R A (q) = v R(q 2 /k 2 ) with R(0) = 1 and R(oo) = 
0. A convenient choice is the Litim regulator R(x) = 
(1 — x)0(l — x), see Ref. [44j . All correlation func- 
tions then depend on the cutoff A. Denoting the flow- 
ing irreducible polarization by IIa(Q), the true irre- 
ducible polarization of our model is recovered in the limit 
limA-»o IIa(<9) = II* (Q). An exact hierarchy of FRG 
flow equations for the one-line irreducible vertices of our 
model can then be obtained by differentiating the corre- 
sponding generating functional T a[(c) , (c) , (<f>)] with re- 
spect to A and expanding T\ in powers of the expec- 
tation values of the fields. A slight complication arises 
from the fact that even in the absence of external sources 
the bosonic field </>q has a finite expectation value </>q , so 
that for finite external sources 

(<t> Q ) = <f>° Q + 6<P Q . (D3) 

In our model, the vacuum expectation value (j>Q in the 
absence of sources is related to the exact density p = 
J k (ckCk) via the Poisson equation^!, 

4>q = 6q,o4> , 4> = -ifop, (D4) 

where Sq^q = f3V8 a .o5 q .Q. The one-line irreducible ver- 
tices can then be defined by expanding the generating 
functional Ta [(c), (c), (</>)] in powers of the expectation 
values %p = (c) , ip = (c) , and 5<fi = (<fi) — <fi° , 



r A [vi,V^ 



yy 1 

^ (m\) 2 n\ 

n=0 m=0 v ; 
i(2m,n) / tsI 



K' 



KL 



Ki 



I<„ 



x T^ m ' n '(K' 1 ,...,K' m ;K 1 ,...,K m ;Q 1 ,. 



<>K[+...+K! m ,K- L +...+K m +Q 1 + ...+Q„ 
i)i>K[ ■ ■ ■ ^K^lpK! ■ ■ ■ 1pK m S(/)Q 1 ■ ■ ■ 5(f)Q n . (D5) 



Note that the vertices T A m ' n> implicitly depend on the vacuum expectation value <j>. Following Ref. [37[ , it is convenient 
to include to contribution (2(3V fo)^ 1 (5(j)a) 2 arising from the fluctuation of the zero mode in the Gaussian part of the 

bosonic action (|3.2p into the definition of the irreducible vertex T^' 2 \—Q, Q) with two external bosonic legs, so that 



-(0,2; 



{-Q,Q) = W^q-o/o -1 +n A (Q). 



(D6) 
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The flowing irreducible polarization ILa(Q) then satisfies the exact flow equatio: 



,36,37 



1 



d A n A (Q) = i/ F A (Q')r^\Q',-Q',Q,-Q)+T^ (Q,-Q,0)d 



(3), 



A 9 A 



F A (Q')F A (Q + Q')r { V(-Q, Q + Q>, -Q')Tf'{Q', -Q - Q', Q)- 



(3), 



(D7) 



Q' 



Here for a sharp momentum-transfer cutoff the bosonic 
propagator is 



F A (Q) = e(q > \q\ > A) : 



u 



i + / ? n A (Q)' 

and the corresponding single-scale propagator is 

U 



F A (Q) = ~S(\q\-A) 



l + f q U A {Q) 



(D8) 



(D9) 



Alternatively, if we work with a smooth additive cutoff 
then 



F A (Q) 



and 



1 + f q [U A (Q) + R A (q)}' 
F A (Q) = {-d A R A (q)][F A (Q)} 2 . 



The vertices 



r^(Qi, ■ ■ ■ , Q n ) — F A ' n> (Qi, 



(D10) 



(Dll) 



(D12) 



r 



are the totally symmetrized one-interaction-line irre- 
ducible vertices with n external bosonic legs. A graphical 
representation of Eq. (|D7[) is shown in Fig. [T6j The flow 
of the vacuum expectation value <p A of the bosonic field 
is determined by the exact FRG equation 37 



[/ - ] - + n A (o)]d A A 



F A {Q')Yf{Q'^Q',0), 



(D13) 

which follows from the requirement that the FRG does 
not generate any tadpole vertices with only one external 
bosonic leg. A graphical representation of Eq. (|D13[) is 
shown in Fig. [I7l We may use Eq. (|D13[) to eliminate 
the derivative of the flowing vacuum expectation value 
in the flow equation (|D7|) to obtain the following exact 
FRG flow equation for the irreducible polarization 



d A U A (Q) = \( F A {Q')Tf(Q',-Q',Q-Q)-\ - - fo — lff>(Q, -Q,0) / F a (Q')T^\q' , —Q', 0) 
2 Jq> 2 1 + /o11a((J) Jq> 

- f F A (Q')F A (Q + Q')T^(-Q,Q + Q',-Q')r^(Q',-Q-Q\Q). (D14) 

Jq 1 

i 



As shown in Ref. [36| (see also the appendix of the 
first work cited in Ref. [37| for a formal proof), in the 
momentum-transfer cutoff scheme the vertices at the ini- 
tial scale A = Ao = qo satisfy non-trivial initial condi- 
tions. The requirement that the vertex T^J 1 '' with one 
bosonic leg vanishes at A = Ao implies that the fermionic 
self-energy £ Ao (A') = T a 2 ^°\k,K) at the initial scale is 
given by the self-consistent Hartree approximation, 

S A = foPo, (D15) 

where the initial density po satisfies the Hartree self- 
consistency condition (|3 . 1 5[) . The initial conditions for 

the purely bosonic vertices r A 3 Q ' ra ' 1 = arc for n = 2, 
ri° ' 2) (-Q,Q) = (PVr^Q.ofo 1 -L ( S 2) (Q,-Q), (D16) 



and for n > 2, 

ri° n) (Qi, . . • , Q n ) = i n (n - 1)! 4 n) H?i> •• • > ~Qn), 

(D17) 

where (Qi, ■ ■ ■ ,Q n ) are the symmetrized closed 
fermion loops with n external legs defined in Eq. (|A2|) . A 
graphical representation of Eq. (|D17|) is shown in Fig. [5] 
By definition, the symmetrized two-loop is (up to a minus 
sign) given by the non-interacting polarization no(Q), 

Lf{-Q,Q) = f G (K)Go(K + Q) = -U (Q). (D18) 

JK 

However, in contrast to Eq. (|2.ip . the fermionic Green 
functions in Eq. (|D18|) are self-consistent Hartree Green 
functions as defined in Eq. (|3.16|) . Finally, in the 
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OX) 



FIG. 16: (Color online) Graphical representation of the exact 
FRG flow equation (|D7[1 for the irreducible polarization in the 
momentum-transfer cutoff scheme. The shaded circles rep- 
resent the one-interaction-line irreducible vertices, the thick 
wavy lines denote the exact cutoff-dependent boson propa- 
gator (effective interaction) defined in Eq. (|D8j) . the small 
crossed circle is the flowing vacuum expectation value of the 
bosonic field <f>, and the small black dot denotes the derivative 
with respect to the flow parameter A. The slashed wavy lines 
represent the single-scale propagator given in Eqs. (|D9ID11[) . 




in) 

mension, then in D dimensions the vertices T A have 
scaling dimensions — (D + z<p)(n/2 — 1), where is the 
dynamic exponent of the bosonic field mediating the for- 
ward scattering interaction. In one dimension z$ = 1 
due to the linear dispersion of the ZS mode, so that in 
D = 1 the vertex r' 3 ' is irrelevant with scaling dimen- 
sion — 1, while r*- 4 ' is irrelevant with scaling dimension 
—2. Because the renormalization group flow of irrele- 
vant couplings is usually not important, it is reasonable 
to truncate the infinite hierarchy of flow equations by 
approximating the vertices and in Eq. (|D7[) by 
their initial values at A = An, 



ri n) (Q! 



) Qn) 



In) 



i n {n-l)\Lf\- 



0, (D20) 



where the symmetrized closed fermion loops are defined 
in Eq. (|A2[) . In particular, for the vertices appearing in 
Eq. (|D14p we substitute 



-^-2iLf{-Q 1 ,-Q 2 ,Q 1 + Q 2 ), 

^a\Qi, —Qi, Q-2, —Q2) 
^6L^\-Q 1 ,Q 1 ,-Q 2 ,Q 2 ). 



(D21) 
(D22) 



FIG. 17: (Color online) Graphical representation of the exact 
FRG flow equation (|D7|I for the vacuum expectation value 
of the bosonic Hubbard-Stratonovich field. The symbols are 
explained in the caption of Fig. 1161 



momentum-transfer cutoff scheme the initial value of the 

C2 1") 

three-legged vertex T K A ' (K'\ K;Q) is 



Tf^iK'-K^Q) 



(D19) 



which follows from Eq. (|3. 3|) . All other vertices vanish at 
the initial scale A = Ao. 

The above RG equation (|D14[) is exact but not 
closed, and should be augmented by FRG flow equa- 
tions for the vertices and , which in turn involve 
higher order bosonic vertices. Only for linearized en- 
ergy dispersion the closed loop theorem guarantees that 

r^CQi, • • • , Qn) =0, so that we recover the well-known 
result that the RPA is exact for the Tomonaga-Luttingcr 
model. To motivate a sensible truncation procedure for 
quadratic energy dispersion, we note that the vertices 
with n > 3 are irrelevant in the renormalization 
group sense: If we assign to the momentum-independent 
part of the interaction /o = / g= o a vanishing scaling di- 



Then Eq. (|D14[) becomes a closed integro-differential 
equation for U\(Q), the solution of which gives for A — > 
a non-perturbative estimate for the irreducible polariza- 
tion. We can easily recover from (|D14[) the perturbative 
expansion (|4. 1|) of the irreducible polarization II* (Q) in 
powers of the RPA interaction if we replace the flowing 
polarization on the right-hand side of Eq. (|D14[) by the 
non-interacting one, Ha(Q) — > Hq(Q), which amounts to 
replacing the flowing effective interaction in F\(Q) by the 
RPA interaction defined in Eq. (|3.23p . With these sub- 
stitutions, it is easy to integrate both sides of Eq. (|D14|) 
over the flow parameter A and recover Eq. (|4.ip . Al- 
ternatively, we may truncate Eq. (|D14[) by replacing the 
flowing polarization on the right-hand side by its limit- 
ing value without cutoff, II* (Q) = limA^onA(Q)- Then 
we arrive that the one-loop self-consistency equation for 

II* (Q) given in Eqs. (|4.13H4.17j) . 

Note that even with the truncation (|D21ID22p the 
FRG flow equation flXUf is non-perturbative, because 
the renormalization of the polarization is self-consistently 
taken into account in the bosonic loop integrations. It 
should be interesting to analyze Eq. (|D14[) numerically 
and try to extract the spectral line-shape. Possibly, one 
can check in this way whether the resummation proce- 
dure proposed by Pustilnik et alA is justified also for 
non-integrable models. 
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